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1 Change of Numeraire 

1.1 Stochastic Financial Modelling 

Recall the following set-up: 

• A market model is a tuple 

M = (n, F, P, (F t )t>o, (Sj\ . . . , S t N ) t > 

where (Cl,J-,¥) is a probability space, (J~t)t a filtration satisfying the usaul conditions, 
and St = (Sf, . . . S^) an (N + l)-dimensional adapted cadlag semimartingale. 

• We will often assume a finite horizon [0,T], e.g. to price European options. 

• We also make the usual assumptions about the market: 

- No transaction costs 

- Continuous trading 

- Liquid markets for every security 

- Short sales allowed 

- Perfect divisibility of assets 

To get results, we will usually specialize: We will generally assume that 0, comes with a 
.ftT-dimensional Brownian motion Wt = (Wj 1 , . . . , W^) which generates the filtration Tt (aug- 
mented). We say that we have K sources of noise. Further, we assume that the asset price 
process is given by an Ito diffusion: 

dS t = M (t, St) dt + a(t, S t ) dW t 

which is shorthand for 



d 



S°\ ffi°(t,S t )\ fo-oi(t,S t ) ... a 0K (t,S t )\ 

dt + 



?N / \ ,,N 



^(t,s t )J 



\<7Nl(t,St) ... 0~Nx(t, St) J 



(dW}\ 
\dW t K ) 
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Under these conditions, the asset price process is (strong) Markov. 

Generally, we make another assumption on S® : we assume that it is the money market 
account process ("riskless" bank account process, which has dynamics 

clSf = rS Q t dt Sg = 1 

A numeraire is a price process N t which has N t > a.s. Think of a numeraire as a unit 
into which other assets are translated. Thus if St is the price of S in money, then St = §^ is 
the price of S in units of N. 

We often choose the numeraire to be the money market account process S®. In that case, 
we write St = for the value of St in terms of the numeraire. Of course, St is just the 
discounted value of S at time t. 

A European contingent claim C is an derivative which, at some future time T has a payoff 
which is a known function of asset prices at time T, i.e. 

Ct = /(St) 

so that Ct is an .^-measurable random variable. The time T is called the maturity or 
exercise time of the claim. 

Our central problem is the pricing and hedging of such derivatives. A European claim can 
be priced by arbitrage methods only if there is a trading strategy which exactly replicates its 
payoff. 

Definition 1.1 A trading strategy/portfolio is a left -continuous (or, more generally, pre- 
dictable) process <pt = (<At\ • • • j which is integrable w.r.t. the semimartingale St- 

</>" is to be thought of as the number of asset S n held in the portfolio at time t. The value 
of the portfolio at time t is 

N 

Vt{<t>) = cj)fSt = Y J <l>tS? 

n=0 

□ 

In the discrete-time setting, we required trading strategies to be self-financing: 

4>t~At ■ S t = 4>f S t or S t A t (p = 

where 4>t is the portfolio held between times t and t+1, and Af(f> = 4>t — <f>t-At- I n continuous- 
time, if we let At — > 0, this looks like an SDE 

S t d<p t = (*) 

However, it would be wrong to use (*) as the self-financing condition in continuous-time, 
because: 

(i) Stochastic integrals are to be interpreted in the Ito sense. 
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(ii) If H t is left-continuous, then the stochastic integral 

fT 

/ H t dX t = lim AX t - X t J 

is a limit (in probability) of left-hand Riemann Stieltjes sums. 

(iii) St((f>t — <j>t-At) looks like a term in a right-hand sum. 

This problem is fixed rather easily: Add and subtract St-At^-t4> from the left-hand side of 
(*) to obtain: 

St-At(<f>t ~ </>t-At) + (St ~ S t - A t)(4>t ~ 4>t-At) = o 
In the continuous-time limit, this looks like 

s t dct> t + d[s,4>]t (**) 

because the first term is a left-hand sum, and the second term looks like a summand in the 
covariation process. 

Thus we would like to define a trading strategy <p to be self-financing if (**) holds. This 
means that the differential of the value process Vt(4>) is easily calculated. By Ito's formula, 

dV t = d(4> t ■ S t ) = 4> t dS t + S t d<f> t + d[S, 4>} t = <h dS t 

and thus ^ 

V t (<t>) = V (4>) + f <P U dS u = Vb(0) + G t (cf>) 
Jo 

where the gains process Gt(4>) is defined to be the obvious stochastic integral. 

Intuitively, this is exactly right: The value of the portfolio at time t equals the initial 
"stake" plus the gain/loss made by playing the stock market game. Nothing else is added (or 
taken away). 

Definition 1.2 A trading strategy 4> is self-financing if and only if Gt(4>) = Jq 4> u dS u , i.e. 
iff d(<j> t ■ S t ) = & • dS t . 

□ 

Actually, the self-financing condition already sneaked into the derivation of the Black- 
Scholes PDE, as given by (e.g.) Hull. There, we started with a portfolio II long one vail option 
C and short n shares S, so that II = C — nS. We then boldly wrote cill = dC — n dS, and 
imposed the condition that II be locally riskless (so that all dWt-terms cancel out, to conclude 
n = But this was dishonest: We should have had dV = dC — n dS — (S dn + d[n, S]), 
because the delta of the option (n) is also changing! However, if we require the hedge to be 
self-financing, then S dn + d[n, S] = 0, and no harm is done. The derivation of the PDE can 
now continue as before. 

Remarks 1.3 In the literature, other conditions are often imposed on trading strategies to 
ensure that they are sufficiently well-behaved. For example, a self-financing trading strategy 
is called tame if Vt(4>) > a.s. It is called admissible if its discounted value is a martingale 
under the EMM. This is important, because even the Black-Scholes model has "doubling" 
strategies, and is not arbitrage-free if arbitrarily large losses can be sustained. However, we 
will ignore thes etechnical points in what follows. 
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□ 

Definition 1.4 A (European) contingent claim C is said to attainable if and only if there 
exists a self-financing strategy (fit such that Ct = Vr(4>) (where T is the exercise date of the 
claim). Then (fi is called a replicating portfolio for C. 

A market model is complete if and only if every contingent claim is attainable. 

□ 

Proposition 1.5 (Numeraire) 

A self-financing portfolio remains self-financing under a change of numeraire. 

This looks totally obvious: After all if we don't add or subtract funds from our portfolio 
when we reckon in units of money, we don't add or subtract funds if we reckon in units of 
barrels of oil either. However, our definition of self-financing is that d{(fit • St) = (fit ■ dSt- Now 
suppose that we reckon in terms of a new numeraire N t . Let St = || be the price of S in units 
of N. To prove that the self-financing condition holds, we must show that d{(fi t • St) = (fit • dSt, 
and this no longer seems so obvious. 

Proof: Let V t = Then by Ito's formula 

^ = w t iV ' +Vtd {w) +i[v ^' 

= |-<is« + *s ( d(i)+*.rf[s,i] 

because (by the self-financing condition) dVt = (fit dSt, so d[V, jj]t = 4>t ■ [S, jj]t- Thus 

= (fitdSt 

□ 

Corollary 1.6 If a contingent claim is attainable in a given numeraire, it is also attainable 
in any other numeraire, and the replicating portfolio is the same. 

□ 

In particular, if the numeraire is the bank account, then 

ft 



Vt(cfi) = V (<t>) + / (fin dS u 
Jo 



Remarks 1.7 A self-financing portfolio cfi = (0°, . . . ,(fi N ) is completely determined by the 
iV of the N + 1 components. Thus if we are given the risky asset components (fi l , . . . , (fi N , 
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the value of the riskless asset component (fP is completely determined by the self-financing 
condition: Take S° to be the numeraire, so that 

n=0 J n =l Jo 

because dS® = — = 1 is constant. Hence 



N \ rt 

n=l U ° 



□ 



1.2 Martingale Pricing 

Definition 1.8 Suppose that N is a numeraire. A measure 
martingale measure (EMM) for numeraire N if an only if 



on (fi, J 7 ) is an equivalent 



(i) 



(ii) S = (j$-)t is a (local) Q-martingale. 
If St is a Q-martingale, Q is called a strong EMM. 

An EMM associated with the money market account is called a riskneutral measure. 



If N is a numeraire, define Vt{4>) = ^7^, and define 



□ 



G 



:(</>)= l\ 

Jo 



dSu 



Note that if Q is an EMM for N, then both V and G are Q~local martingales. Indeed, 
Gt = Jq <fiu dS u is a sum of stochastic integrals w.r.t. a Q-local martingale. 

We require Q to be equivalent to P so that both measures have the same arbitrage strate- 
gies: F(G T > 0) > if and only if Q(G T > 0) > 0. 

Further note that 

• An arbitrage opportunity remains an arbitrage under 

— a change of equivalent measure; 

— a change of numeraire. 

• A replicating portfolio remains a replicating portfolio under 

— a change of equivalent measure; 

— a change of numeraire. 

Theorem 1.9 // an EMM Q exists (for some numeraire N), then there are no arbitrage 
opportunities. 
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Proof: If <f> is a self-financing strategy, then 

= G (<l>) = Eq[G t (<I>)] 



Now because P(Gt > 0) > if and only if Q(Gt > 0) > 0, and because G T > if and only 
if Gt > 0, we cannot have both Gt > and Ep[Gt] > 0. Thus (j) cannot be an arbitrage, i.e. 
there are no arbitrage opportunities. 

□ 

Example 1.10 The most common choice of numeraire is the money market account. Sup- 
pose that 5° is the MMA, with price dynamics 

dS? = r(t,u)S$ dt 

If Q is the EMM associated with S°, then each S™ is a Q- local martingale. Now 

iqn on 



and so 



dS? = S°dS? + rS™ dt = rS? dt + dM? 



where M t n = J * Sf dSf is a Q-local martingale. Conversely, if each dS? = rSf dt + dM? for 
some Q-local martingale M™, then Q is a riskneutral measure. 

□ 

Theorem 1.11 (Riskneutral Valuation) 

Suppose that X is an attainable contingent claim, and that Q is an EMM for numeraire N. 
Then 

X t = EQ[X T \Ft] 



i.e. 



X t = N t E Q 



Xt 



1^ 



Proof: If (f> replicates X, it does so under any numeraire, any EMM. Now by the Law of One 
Price, 

X t = V t {<t>) = E Q [V T (cl>)] = E Q [X T ] 

□ 



1.3 Introduction to Change of Numeraire 

Thus far, we've mainly considered two probability measures, the "real world" measure P, and 
the equivalent martingale measure Q for the money market account numeraire. We've seen, 
however, that it is possible to introduce an EMM for different numeraires, and to use these 
for pricing. We now show that a change of numeraire is a technique which often simplifies a 
pricing problem - it is analogous to a reduction in dimension. 
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Consider an interest rate derivative X, and let A t be the bank account. If Q is the EMM 
for A t , then ^ is a Q-martingale (assuming, of course, that X is attainable). Thus 



An = E f 



e -/ T r(s,u;) ds x ^ 



where r is the short rate, so that A t = e-fo r ( s > w ) ds . I n order to compute this, we would have 
to know the joint density of X(T), A(T) under Q — it would not be observable, because only 
P-densities can be observed. The computation of the expectation would involve a double 
integral. 

The reason we haven't noticed this problem before is that we generally assumed that 
interest rates are constant, which simplifies matters considerably. If we assume that the 
payoff X(T) and the short rate are independent under Q, then we would still have some 
simplification, namely 



X n = E f 



- Jo r{s,w) ds 



= p(0,T)E Q [X(T)\ 

where p(t, T) is the time t-price of a zero coupon bond with face value 1 and maturity T, 
i.e. an interest rate derivative with payoff 1 at expiry, in all states of the world. The above 
expression is obviously much simpler: 

• It only involves a single integral, and needs only the Q-density of X{T). 

• p(0,T) is observable (either directly, or by bootstrapping a yield curve from observable 
coupon bond prices). 

Generally, of course, X(T),A(T) are not independent under Q. Even if they were independent 
under P, they would nevertheless probably not be independent under Q — under Q, the drifts 
of all assets are the same, namely the short rate. Thus X t has the same drift as A t , implying 
some correlation. 

1.4 Mechanics of Changes of Numeraire 

As usual, we work with a market model (O, T , P, {J~t)t-, ■ ■ ■ , S 1 /^)*). Recall: 

• A numeraire is a traded asset (posibly a portfolio of assets) with a strictly positive price 
process. 

• Self-financing portfolios remain self-financing under a change of numeraire. 

• Replicating portfolios remain replicating portfolios. 
This is also a good time to recall: 

Theorem 1.12 (Bayes' Theorem) 

Suppose that (fi, T, P) is a probability space equipped with a filtration T n , and that Q « P. 
Let £ = dQ/dP and likelihood process = Ep[£|J^]. Then 

(a) For any random variable Z (integrable w.r.t. P and Q) we have 
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(b) If Q « P, i/jen a stochastic process X t is a martingale under Q if and only if £ t X t is a 
martingale under P. 

The proof is an exercise: 

Exercises 1.13 Suppose that (Tl,J-,F) is a probability space equipped with a filtration J- n , 
and that Q « P. Let £ = dQ/dF and define ^ = Ep[£|.Fi], where Ep refers to expectation 
w.r.t. the measure P. 

(a) Show that for any random variable Z (integrable w.r.t. P and Q) we have 

(b) Show that if Q « P, then a stochastic process Xt is a martingale under Q if and only if 
$,tXt is a martingale under P. 

[Hint: (a) I'll give you the proof. You justify every step: Let A G T n . Then 

it^\Z\T t \dF= [ E P [ZE Q [Z\F t }\F t ] dF 

A J A 



J A 

[ Z dQ 

J A 

[ ZidF 
J A 



I A 

(b) Use (a). ] 

□ 

Now the bank account Sf = At is just a special numeraire — one whose dynamics have zero 
volatility: dA t = r(t,u)A t dt. Let Q be the EMM for A t . Then each -fc- is a Q-martingale, 
i.e. Q "martingalizes" the ratios 

Suppose that A(t) is another numeraire, with EMM Q. Q "martingalizes" the ratios -4-. 

At 

Given that A{t) is a (combination of) traded assets, we expect to be a Q-martingale as 

well. If X is an attainable claim, then ^ is a Q-martingale, and ^ is a Q-martingale. 

What does Q look like? Since Q,Q are both equivalent to P, they are equivalent to each 
other, and thus the Radon-Nikodym derivative 

m = d 4 

dQ 

exists. We don't know yet what it is, though, because we don't know Q. Nevertheless, it 
exists, so we may define the likelihood process 



L(t) = E Q [L(T)\F t ] 
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We have, by Bayes' Theorem, 



io 



X(T) 
A{T) 



= L(0)- l Eq 



X{T) 
[A(T) 



L{T) 



so that 



because L(0) = 1. But 



A 



Er 



X(T) A 



A n 



E 



X(T) 
A(T) 



L(T) 



This suggests that we turn every thing around and define 



L(T) 



A{T)/A(T) 



and then define Q by 



L(T) 



Then L(t) = , as y°u can easily check. 

In general, we may use for A(t) absolutely any process with the property that ^ is a 
strictly positive Q-martingale. 

Theorem 1.14 (Martingale Measure Pricing) 

Suppose that A{t) is process with the property that ^ is a strictly positive Q-martingale. 
Define 

A t /A t 



-(f i> s a Q-martingale, then ^ is a Q-martingale. In particular, if X is an attainable 
contingent claim, then 

r X(T), 



L(t) 



Mt 



Aq/Aq 



L(T) 



X t = A t E f 



A(T) 



□ 



In fact, we can generalize even more: 

Theorem 1.15 Suppose that ai(t) , a 2 {t) are numeraires, and that Qi , Q2 are their associated 
EMM's. Then for any random variable X we have 



ai(t)E Ql 



X 
ai(T) 



a 2 (t)E Q2 



X 
a 2 (T) 
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Proof: Define the likelihood process L\{t) = 2lfo)/MQ) > anc ^ define L,2(t) similarly. Then by 
Bayes' Theorem 



Qi(t)E (Ql 



and the same goes for «2- 



X 



a 1 (t)L 1 (t)- 1 l 
A(t)E< 



X 



«i(T) 



Li(T)\F t 



A{T) 



□ 



Let's investigate how asset price dynamics change when we move from Q- world to Q- 
world. Assuming that asset prices are Ito diffusions, we have Q-dynamics 

dS t = D[S t ]r t dt + D[S t ]a t dW t 

where W is a (if-dimensional) Q-Brownian motion. The Radon-Nikodym derivative (pro- 
cess) which effects the change from Q to Q is 



L{t) 



Mt)/A(t) 

i(o)M(o) 



Using the fact that A t /A t is a Q-martingale, we see that 

dA t = r t A t dt + a t A t dW t 

By Ito's formula, 

A 



dL t 



An 



[r t A t dt + a t A t dW^j - ^ [r t A t dt) 



Thus 



dL t = L t a t dW t 

confirming that L(t) is a Q-martingale, as we already knew. Solving the SDE, we obtain 

L ( t ) = e /> s dw s -y*\\a s \\ 2 ds 

and thus 



Suppose we change the numeraire from the MMA At to At. Then the 
EMM Q associated with At is obtained from the EMM Q associated 
with At by a Girsanov transformation whose kernel is the volatility a of 
the new numeraire At. 



By Girsanov's Theorem 



t 

W t = W t - I a s ds 
h 



12 



Change of Numeraire 



is a (K-dimensional) Q-Brownian motion, and thus the asset dynamics under Q are given by 

/ r t + a] ■ a t \ 



dS t = D[S t ] 



dt + D[S t ]a dW t 



\r t + a? ■ a t J 



i.e. 

dS? = {r t + cj t n • o)S? dt + a?S? dW t 
where a n is the n th row of the volatility matrix a. In particular, the "discounted" asset ratios 
St = ~A have dynamics 

dS^Sno-Z-aJdW* 

as you can easily verify by applying Ito's formula. Hence the S^ are Q-martingales, and the 
volatility of each "discounted" asset is reduced by the volatility of the numeraire. 

Remarks 1.16 Consider a simple Black-Scholes model, where the risky asset prices are given 
by a geometric Brownian motions, driven by a single source of noise. The market price of 
risk of S n in Q-world is 

r + a n a — r 

= a 

a n 

i.e. all assets have the same market price of risk, namely the volatility of the numeraire. This 
is also true in the multidimensional case, where the market price of risk is a vector. 
The bank account has zero volatility, and thus the MPR in Q-world is zero. 

□ 



1.5 A General Option pricing Formula 



Consider a call C on a security S with strike K and maturity T. Let Qs be the EMM 
associated with numeraire S, and let Q be the T-forward measure (i.e. the EMM associated 
with the zero coupon bond p(t,T) maturing at T). 



Theorem 1.17 
Proof: We have 



C = S q s (S T >K)- KP(0, 



T 



(S T > K) 



C = p(0, T)E Q T[max{S T - K, 0}] 
= p(0,T)E QT [S T -K;S T >K] 

= p(0, T)E QT [S T ; S t >K\- Kp(0, T)Q T (S T > K) 



But we have 



X 



ai(T) 



\Tt 



X 



a 2 (T) 



for general numeraires and their associated EMM's. Using this with a\(t) = p(t,T) and 
&2(t) = St, we obtain 



p(0, T)E QT [S T ; S t >K}= p(0, T)E q t 



SqEq £ 



StI{s t >k} 
P(T,T) 

StI{s t >k} 



St 

= S <Qs(S T > K) 
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□ 



1.6 Applications 

Example 1.18 Forward Measures 

Consider again the situation at the beginning of this chapter: We consider a contingent claim 
X with expiry T. Under riskneutral valuation, its value is 

X = E Q \ e ~So r M ds X(T) 

where r is the short rate. We bemoaned the fact that this would necesssitate us knowing the 
joint density of At and Xt- If only, we said, At and Xt were independent, we would get the 
much simpler 



Xn 



- Jo r(s,w) ds 



®q[x(t)] 



= p(0,T)E Q [X(T)\ 

where p(t, T) is the time t-price of a zero coupon bond with face value 1 and maturity T. If 
only. . . 

Now let's see what happens if we change the numeraire to p(t,T), and let Qt be the 
corresponding EMM. In that case, the pricing formula becomes 



Xn 



Xt 



[p(T,T) 



p(0,T) 

and noting that P(T, T) = 1, we have 

X o = p (0,T)E Qt [X t ] 

This is the simple form that we sought, but it's correct under Qt, and not under Q. 

The measure Qt is called the T-forward measure. Note that if interest rates are de- 
terministic, then Q and Qt coincide, because then p(t,T) = e~ $t Ts ds 7 so that each ratio 
S?/p(t,T) = A T (S?/A t ) = const, x S?/A t is already a Q-martingale. 

However, when interest rates are stochastic, Q and Qt are quite different. We shall see 
later that futures prices are Q-martingales, whereas forward prices are QT-martingales. Thus 
forward prices and futures prices coincide if interest rates are deterministic. 

□ 

Example 1.19 Exchange Options 

Consider an exchange option which gives the right, but not the obligation, to exchange asset 
S l for asset S 2 at time T. This is a contingent claim X with payoff 

X T = m&x{S T - S T ,0} 

Using riskneutral valuation (i.e. MMA as numeraire), its value is therefore 



e-^ rtdt max{S T - S T ,0} 



To compute this, we have to know the joint distributions of At,S t ,S t under Q, yielding a 
triple integral. 
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It is computationally simpler to change the numeraire: Let A t = S^, and let Q be the 
associated EMM. The contingent cliam is then priced as follows: 



1-2 - S? 



max{Sj. — S^, 0} 



max{^ - 1,0} 



(where Sf = ^). This looks like a call on S 1 with strike K = 1, and we only have to know 

the distribution of Sj. under Q. 

To price this option, we have to assume some form of asset dynamics. Suppose these are 
given by one-dimensional Ito diffusions, i.e. suppose we have P-dynamics 

dS] = inS} alt + cti dW^(t) 
dSl = fi 2 S? dt + a 2 dWp(t) 

where Wp,Wp are correlated P-Brownian motions, with correlation p t . To get lognormality 
we also assume that d\{t), a 2 (t) and p(t) are deterministic. Interest rates, however, can be 
stochastic. 

We can write this as 



s 2 t ) U$ 



where Wl , Wp are independent P-Brownian motions and 



a 



0~11 CT12 
C21 C"22 



Then we must have 



o x = yj a\ x + a 2 



12 



a 2 = \/crfi + <? 22 and 

0~llCr 2 i + CT12CT22 



p 



^11+^12^21+4 



22 



So given c7i,<72 and/? we can solve for the matrix a (though not uniquely). 

Note that the correlation s a function of the volatility matrix. When we change from 
P-world to Q-world, the volatility matrix is unchanged, and thus also the correlation. Thus 
under Q, the asset dynamics are 



dSl = r t S] dt + a x dW^(t) 
dS? = r t S? dt + a 2 dW%(t) 



where WX, Wq are Q-Brownian motions, with correlation p. This can also be written as 



d ( S Jf ) =D[S t ]r t dt + D[S t ]a t 
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where Wq,Wq are independent Q-Brownian motions. Now when we change the numeraire 
from the MMA to S 1 , and the measure from Q to Q, we get 

fWl (t) 

dSt = S?(a 2 i - an, (7 2 2 - 0-12) • d ® 

\W 2 (t) 

= S 2 a 2 dW^t) 
where Wq is a one-dimensional Q-Brownian motion and 

0-2 = V (021 - (Tn) 2 + (o"22 - 0-12) 2 = \J (^l) 2 + (02) 2 - 2po\o 2 

Now since the o\ (t) , a 2 (t) and p(t) are assumed to be deterministic, so is a 2 (t). It follows 
that is lognormal under Q: Indeed 

rT ~n\ jw. u\ 1 r T s. 



so that 



ln(S 2 /S 2 )~iv(-^ a 2 (t) 2 dt, 1 -J^ a 2 (t) 2 dt 
Using the properties of lognormality, we see that 



where 



max{5^ -1,0} = E 6 [S%\N(d{) - IN fa) 



di = ln(^) + i/ T a 2 (t) 2 dt 



dt 



d 2 = d x - ^ a 2 {t) 2 dt 

and thus, using the fact that Eq[5j>] = S 2 , we have 

X = SiNfa) - S^N(d 2 ) 
If we further assume that a±,a 2 and p are constants, we obtain 

X = 5 2 iV(di) - S*N(d 2 ) 



where d\ 



ln(f ) + l(a 2 + a 2 - 2pa 1 a 2 )T 



^a 2 + a 2 - 2pa 1 a 2 )T 
d 2 = di - \Jv\ + dl- 2pa\o 2 )T 

where we used the fact that 

a 2 = v 7 {°i) 2 + (& 2 ) 2 - 2pa\a 2 

□ 
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Note that, in the above example, ^ J Q T <5"2(i) 2 dt is just the average of the squared volatility 
of S 2 , so that J Q T a 2 {t) 2 dt = o- 2 vcragc T. 

Example 1.20 A convertible bond is a bond (issued by a corporation) that can be converted 
to equity at certain times, using a predetermined exchange ratio. For example, a bond with 
a conversion ratio of 10 allows its holder to convert a par $1000 bond to 10 shares of the 
common stock, at some future date T (which may, or may not, e the maturity of the bond). 
Thus 

^ - n i n ■ n - n Bond price 
Convertible bond = regular bond + 10 calls with strike — 

Here the bond price is the future bond price at the convertion date T, which is unknown if T 
is not also the maturity of the bond. The strike varies because the bond price varies. 

Thus a convertible bond will be difficult to price if interest rates are stochastic, and 
practitioners will generally use an estimate of the future bond price. For example, if the 
bond's coupon rate is approximately equal to the prevailing interest rate, then setting the 
future bond price equal to its par value will not be a bad approximation. The fact that the 
corporation typically issues new shares to honour the conversion also means that a dilution 
effect must be taken into account. 

We make the following simplifying assumptions: 

• The bond is a zero coupon bond, with face = 1, maturity T±; 

• The underlying stock pays no dividends; 

• At a fixed date Tq < Ti, the bond can be converted to cb shares of stock. 

Our aim is to price this convertible bond at some time t <Tq. 
The payoff of this bond is 

X To =mzx{cS To ,p(To,T 1 )} 

where p(t,T\) is a face = 1 zero coupon bond of the same risk class as the convertible bond. 

In order to simplify this expression, we should use either S or p(t,T\) as numeraire. In 
the first case, we obtain 



j>(T ,Ti) 

maxjc, — }\F t 

b T 



max{ l}\F t 



X t = StE(Q s 
whereas the second choice yields 

X t = p(t,T 1 )E QTl 

L b 

Let's consider the second possibility: Define St = ^fjyj- Then 

X t = cp^T^Eqt, [max{5 To - l/c,0}|^ t ] +p(t,Ti) 
Now we can observe the market value of p(t, Ti), so we need only calculate 

cp(t,Ti)EQ Tl max{S To - 1/c, 0}|^i 
which is just the price of c calls with strike K 
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Now in order to price the option, we need to know something about the distribution of 
St - This requires the specification of both a stock price model and an interest rate model. 
For the stock price model, we can take the usual GBM. As for interest rate models, this is 
the subject we will tackle next. We will continue this example in the Exercises, where you 
will be asked to price a convertible zero-coupon bond in the Ho-Lee model. 

□ 

2 Modelling Fixed Income: Introduction 

2.1 Classification of Interest Rate Models 

We will examine several approaches for the modelling of interest rates: short rate modelling, 
whole yield curve modelling and market models. The purpose of this section is to introduce 
basic concepts and notation. Amongst the immediately obvious quantities that we may model 
are 

• bond prices 

• the short rate 

• forward rates (discretely or continuously compounded) 

• the entire yield curve 

Of course, a model of bond prices will have the yield curve as an output, etc. These approaches 
are no independent. 

Short rate models: These model just one variable, the short rate, which is an idealized 
quantity that represents the instantaneous interest rate at any time. Usually a diffusion 
model, and thus Markov. We specify dynamics, e.g. 

dr t = k(9 — r t ) dt + a dWt k, 6, a const. 

is the Vasicek model, and 

dr t = n(0(t) — r t ) dt + a(t) dWt k const., 6, a deterministic 

is the Hull-White (extended Vasicek) model. 

• Can be one-factor or multi-factor 

• Affine term structure models have a particularly simple form, allowing for closed form 
solutions for bond option prices, Eurodollar futures, etc. More later. . . 

• Multi-factor models: principal component analysis shows that 80-90% of the variance 
of the term structure is explained by parallel shifts of the yield curve, 5-10% by a twist 
(long term and short term rates move in opposite directions, pivoting about a point), 
and 1-2% by a butterfly (long and short term rates move in the same direction, with 
mid-term rates moving in the opposite direction) . 
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Whole yield curve models: These model the entire term structure of rates, eg. the entire 
forward rate curve. Examples are 

• Heath-Jarrow-Morton models 

• Market models 

Interest rate models are often categorized into Equilibrium models and No-arbitrage mod- 
els. Equilibrium models attempt to derive, e.g., short rate dynamics from macroeconomic 
considerations, starting from a representative investor (e.g. Cox-Ingersoll-Ross, Vasicek, 
Merton models). These models often have the nice property of being time-homogeneous, 
but usually are unable to fit observed prices exactly. No-arbitrage models attempt to fit a 
model exactly to observed prices and volatilities -zero coupon bonds, caplets, swaptions. (e.g. 
Ho-Lee, Hull- White models). 

Both terms are misnomers: Some equilibrium models are not arbitrage-free, and thus 
not in equilibrium. Some no-arbitrage models permit negative interest rates, thus allowing 
"mattress arbitrage" (borrow from the bank when rates go negative, put under mattress). 

2.2 Bond Market Basics 

One of the basic instruments that we shall be concerned with is the following: 

Definition 2.1 A T-bond is a zero coupon bond with face value 1.00 and maturity T. Its 
value at time t < T is denoted by p(t, T). 

□ 

These are also called discount bonds. 

We work in a probability space (£l,J-,F) equipped with a filtration which satisfies the 
usual conditions. We usually require that: 

• p(t, T) is a continuous semimartingale for each T. 

• < p(t, T) < 1 a.s. (This fails in, e.g. Gaussian short rate models) 

• There is a frictionless market for T-bonds of every maturity T > 0. 

• For every fixed T > 0, {p(t, T) : < t < T} is an optional process with P(T, T) = 1. 

• For every fixed t, p(t,T) is (P-a.s.) differentiable in the second variable T. 

,. ^ dp(t,T) 

• No default risk. 

Note that, for fixed t, the set {p(t,T) : T > 0} is just the term structure of zero coupon 
bond prices, which is typically a smooth decreasing function (of T). On the other hand, for 
fixed T, the set {p(t, T) : t < T} is the price process of the security p(t, T), which is typically 
very ragged (i.e. of unbounded variation). 

Note that there are, in our model, infinitely many securities, namely one p(t, T) for each 
maturity T. 

We briefly recall the definitions of the various types of rate: 
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• Let t < S <T. Consider the following strategy: 



(i) At time t, short an S'-bond, and use the proceeds to buy ^j^ -many T-bonds. 
Net cashflow at time t is zero. 

(ii) At time S, pay $1.00 to redeem the 5-bond. 

(iii) At time T, receive ^'^j from maturing T-bonds. 

Thus at time T, we can, with no initial cash outlay, ensure that a deposit of $1.00 at 
time S leads to a payoff of ^[pfy at time T. This implies that we can lock in an interest 
rate R(t; S, T) for the future period [S, T]: 

p R(T-S) = gftvg) 

~p(t,T) 

^R(t;S,T) = - lnP ^ T l- l ^ S) 

1 — o 

This is the forward rate (continuously compounded) for the period [S, T] at time t. 

• The equivalent simple forward rate (the LIBOR forward rate) for [S, T] contracted at 
time t is similarly defined by 

p(t,S) 



1 + L(T — S) 

=>L(t;S,T) 



p(t,T) 
p(t,T)-p(t,S) 
p(t,T)(T-S) 



• The continuously and simple spot rates at time t for time T are R(t; t, T) and L(t; t, T) 
respectively. 

• The instantaneous forward rate at time t for time T is the interest rate that can be 
locked in for an infinitesimal interval [t, T + dT]. It is given by 

f(t, T) = R(t; T, T + AT) 

_ Um lnp(t,T + AT) -]np(t,T) 
~ A™o AT 

dlnp(t,T) 
~ df 

• The short rate is the instantaneous spot rate, and is defined by 

r(t) = f(t,t) 

• Given a tenor structure 

< t < T < Ti < T 2 < • • ■ < T N 

we can find a forward swap rate St = S(t; Tq,T±, . . . , T/v), the unique fixed rate, at time 
t, for which a fixed-for-floating forward swap, starting at time To, will have zero value. 
We clearly require, with t 3 - = Tj — Tj-i, that 

N N 

s t T jP (t, Tj) = m Tj-utyTjpit, Tj) 

3=1 3=1 
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and thus 

ZUmTj^T^rjp^Tj) 

St = — 



But 

N N 

L(t; Tj- U Tj)Tjp{t, Tj) = ~\p(t, T j) ~ P(t, Tj-i)} = p(t, T ) - p(t, T N ) 

3=1 3=1 

and hence 

_ p(t,T ) - p(t,T N ) 
St - 



The denominator J2jLi T jP(t,Tj) is sometimes referred to as the value of a basis point. 

Remarks 2.2 1. The assumption that there are traded zero coupon bonds of every maturity 
is clearly false. Nevertheless, a large number of implied zero coupon bond prices can usually 
be obtained by bootstrapping the yield curve. 

2. The instantaneous rates (forward- and short-) are theoretical entities, and not directly 
observable in the market. One of the shortcomings of short rate and HJM models is that 
they model these non-existent entities. Market models such as the BGM- and Jamshidian 
models, however, are concerned with the modelling of quoted market rates. 

□ 

The following lemma shows how bond prices are related to forward rates: 
Lemma 2.3 

p(t,T)=p(t,S)e-fsf( t > u ) du 
Proof: In p(t, T) = In p(t, S) + /J ^ du. 

As usual, we denote that money market account (MMA) process At by 

A t = ef° r(M) du 

where r(t) is the short rate. 

Example 2.4 No model that allows only parallel shifts of the yield curve is arbitrage-free. 
Proof: Suppose it is certain that f(l,T) = f(0,T) + e for all T > 1, where e is a random 
variable. Now choose times 1 < 7\ < T 2 < T 3 . At t = 1, 

v (l T ) - e - £ f{l,v) du _ - f(0,u)+e du _ gCjVO -e(T-l) 

Pl ' j ~P(0,1) 

Now suppose that we hold Xj Tj-bonds (i = 1,2,3). We construct an arbitrage, a static 
portfolio satisfying 



□ 
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(i) EL^(0,T J )=0 

(ii) Eti^P(l,^) >Oa.s. 

At time 1 the value of the portfolio is 

3 



z=i 

3 

= £ 

i=l 
3 

-E 



i=l 
3 

Xz ~ 

i=l 
3 



p(Q,ri) r - em -i) 

p(0,l) 



g g(0 1 ji) e - e (r < -r 2 ) e - e (r 3 -i) 



, : 1*0,1) 

e -e(2b-l) 

9<E) iH<ur 



where 



3 

=^^(0,^)6-^-^) 
i=i 

We shall ensure that Vi(e) > whenever e ^ 0. First note that g(0) = 0, because 
Z)i=i x iP{^i Ti) = 0. Further, V\{e) and 51(e) always have the same sign, so to ensure V\{e) > 0, 
it suffices to ensure that g(e) > 0. 

Now g is a C 2 -function (twice differentiable) , and we require that (i) 51(0) = 0, (ii) g(e) > 
whenever e / 0. It follows that g'(0) = 0, thus that 

3 

g'{Q) = Y J X i {T 2 -T i )p{Q,T i ) = Q 

i=l 

and thus that 

3 

i=i 

Next, if we ensure g"(e) > 0, then, combined with g(0) = <?'(0) = 0, we se that g(e) > 
for all e / 0. Now 

g"(e) = J>(T 2 - TOyO.TOe-^- 2 *) 
i=i 

and thus g"(e) > for all e if xi, X3 > (and at least one is > 0). 

Now take x 2 < 0. Since Ylf=i x iP(®i ^i) = 0) we see that at least one of xi,X3 must be 
> 0. Since Ei=i x «(^2 — Ti)p(0, Tj)p(0, Tj) = 0, we see that xi,X3 have the same sign, i.e. 
both are > 0. Then g"(e) > for all e / 0, and hence also g(e) > 0. 

It follows that any portfolio (xi,x 2 ,xs) satisfying 

W Eti^p(o,r,) = o 

(ii) E 3 i =i x i( T 2-Ti)p(0,T i ) = 

(iii) x 2 < 
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is an arbitrage. 

□ 



Example 2.5 Define the long rate l(t) by 

l(t) = lim R(t, T) 

where R(t,T) is the c.c. spot rate, i.e. p(t,T) = e~ R ^' T \T — t). Though l(t) is not directly 
obtainable form traded securities (because the longest-term securities typically have a life 
of 30 years or so), it can be estimated, and empirical studies suggest that it fluctuates con- 
siderably over time. Most no-arbitrage models have a constant value for l(t), however, and 
indeed 

Theorem:// the term-structure dynamics are arbitrage-free, then l(t) is an increasing func- 
tion a.s. 

Proof: By rescaling time, we may assume that 1(1) < 1(0) with positive probability, to 
obtain a contradiction. For T = 1,2,3,..., construct a portfolio which, at t = invests 
^ — = T( - r 1 +1 ^ into each of the bonds p(t, T), so that the value of then portfolio is Vq = 

Et=i TiT+T) = L Define e = W°) " Z ( 1 ))/ 3 - Now P(°> T ) = e- r (°> T ) T , and r(0,T) -► 1(0) 
as T — > oo, so eventually, we have r(0,T) > 1(0, T) — e, i.e. p(0,T) < e~W)~ £ ) T eventually. 
Similarly, p(l,T) > e~^^ +e ^ T eventually. Suppose these relations hold for all T > Tq. Then 

v = y P(1,T) ^ P (1,T) 

1 2^ 7YT 4- 1 WO T~l 7YT 4- 1 Wfl T) 



- T(T + l)p(0, T) ^ T(T + l)p(0, T) ^ T(T + 1) 

The second term diverges to oo, so that V\ = oo. Now since Vo = Eq[Vi/.Bi], where Q is 
a risk-neutral measure and B is the bank account, we see that Q(V"i = oo) = 0, because 
Vq = 1 < oo. Since the "real-world" measure P is equivalent to Q, we must have P(Vi = 
oo) < P(Z(1) > 1(0)) = as well. 



□ 



2.3 Modelling the Bond Market 

We consider three approaches: 



1. Specify short rate dynamics; 

2. Specify bond price dynamics; 

3. Specify forward rate dynamics; 

Suppose, for example, that we are given the following dynamics: 
1. Short rate dynamics: 

dr(t, uj) = a(t, uj) dt + b(t, uj) dWt 
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2. Bond price dynamics: 



dp(t, T)(lo)= p(t, T) (u) [m(t, T, uj) dt + v(t, T, u) dW t ] 



3. Forward rate dynamics: 



df(t,T)(co) 



a(t, T, u) dt + a(t, T, uj) dW t 



Here Wt is a standard (multidimensional) Brownian motion. 

If we're given one type of dynamics, can we deduce the others? If you think about this for 
a while, you'd expect that bond prices and short rates are deduceable from the forward rates, 
and that forward rates and the short rate are deduceable from the bond prices. A model 
of just the short rate seems to contain too little information to deduce all bond prices and 
forward rates however. 

Before we write down exactly how the various dynamics are related to each other, we need 
a stochastic Fubini Theorem and its corollary. 

Proposition 2.6 (Fubini's Theorem for Stochastic Integrals) 



where (s,cu,S) t— ► <&(s,S,uj) is V x B-measurable (V = predictable a-algebra, B = Borel 
algebra), and 



The proof is omitted, but may be found in Durrett, Chapter 2, Section 11. 

Before we prove a corollary about the differentiation of stochastic integrals, it is convenient 
to gather well-known results about the differentiation of ordinary Lebesgue integrals: 

Proposition 2.7 Assuming sufficient smoothness and regularity, 





□ 




d_ 

dx 



rh(x) rh[x) Q dh 

/ f(x,y)dy= —f(x,y)dy + f(x,h(x))—-f(x,g(x)) 



dx 



dg 




□ 
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Corollary 2.8 (Differentiation under the integral sign) 

r * rt dv(s,T) 



d_ 
df 



f v(s,T)dW s = [ 

J 



dT 



dW, 



Proof: Just like the ordinary proof of differentiation under the integral sign: 





dT 



fv{s,T) 
Jo 



dW, 



d_ 
df 



v(s,0) + 



dv(s, u) 



dT 



du dW s 



+ 



t r T 



d_ 

dT J J 



dv(s, u) 



dT 



du dW« 



--if 

9TJo Jo 

-f 
Jo 



dv(s, u) 
dT 



dW* du 



dv(s,T) 
dT 



dW, 



□ 



Consider now the various dynamics given above, i.e. short rate, bond price and forward 
rate dynamics. Assume that the drifts and variance rates are C 1 in the T-variable, and 
sufficiently regular to allow the interchange of order of integration. Further, assume that 
bond prices are bounded. 

The following theorem records the relationships between the various dynamics: 



Theorem 2.9 (a) If 

then 
where 

(b) U 
then 
where 

(c) If 
then 



dp 
V 



= m dt + v dWt 



df = adt + a dW t 

a(t, T) = v T (t, T) ■ v(t, T) - m T (t, T) 
a(t,T) = -v T (t,T) 

df = adt + a dW t 
dr = adt + b dW t 

a(t) = f T (t,t)+a{t,t) 
b(t) = a(t,t) 

df = adt + a dW t 



dp 
p 



r(t)+A(t,T) + -\\S(t,T)\\ 2 



dt + S(t, T) dW t 
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where 

A(t,T) = -J a(t,s)ds 

S(t,T) = ~J ds 
Here \ \ ■ \ \ is just the usual Euclidean norm. 

Before we begin the proof, note that for each T we have a separate security p(t,T), i.e. 
for each T we have a separate process (p(t, T))t>o- It is to these processes that we apply Ito's 
formula, etc. 

Proof: (1) dlnp = [m — \v 2 } dt + v dWt and thus 

lnp(t,T) = lnp(0,T) + f m(s,T) - 1 v 2 (s,T) dt + f v(s,T) dW s 

Jo 2 j 

so that 

„ ^ d\np(t,T) d\np(0,T) /"* , /"* 

-f(t,T) = ^-L = + J^m T - VT . v ds + J v T dW s 

Taking differentials yields the result. 
(2) 

r(t) = f(t,t) = f(0,t) +/ a{s,t)ds + a(s,t) dW s where 

Jo Jo 

a(s,t) = a(s, s) + J ax(s,u) du 

cr(s,t) = cr(s,s) + J gt{s,u) du 

and thus 

rt rt rt rt rt rt 

r(i) = /(0,t) + / a(s,s)ds + / / ar(s,u)duds + I a(s, s) dW s + / / ctt(s, u) du dW s 
Jo Jo J s Jo Jo Js 

rt rt ru rt rt ru 

= f(0,t)+ / a(s,s)ds + / / ar(u, s) ds du + / a(s,s)dW s + / / <tt(u, s) dW s du 
Jo Jo Jo Jo Jo Jo 

by the stochastic Fubini theorem. Thus 



dr(t) 



a(t,t)+ I a T (s,t)ds + f a T (s,t) dW s dt + a(t,t) dW t 
Jo Jo 

a{t, t) + fr(t, t)} dt + a(t, t) dW t 



as required 

(3) First define Y(t,T) = - / T f(t,s) ds, so that p{t,T) = e Y ^ T \ Now 



f(t,s) = f(0,s) + [ a(u,s)du+ [ a(u,s) dW s 
Jo Jo 
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and hence 



Y(t,T) = - f(0,s)ds- / a(u,s)duds- / o{u, s) dW u ds 
Jt Jt Jo Jt Jo 

T ft fT ft fT 

/(O, s) ds — \ \ a(u, s) ds du — J / a(u, s) ds dW u 
Jo Jt Jo Jt 

T ft 



- / f(0,s) ds+ f f(0,s) ds 
Jo Jo 



ds du 



+ 



t rT rt rt 

/ a(u, s) ds du + / / a(u,s) 

O J u Jo Ju 

t fT rt rt 

u, s) ds dW u + / o-(u, s) ds dW u 

Jo Ju 



-[ I a{r 

JO Ju 



Y(0,T) 



t 



a(u, s) ds du 



o Ju 
t 



t 



a(u, s) ds dW u 



O Ju 



ft ft ft ft ft 

+ / /(O, s) ds + / a(u, s) ds du + / o~(u, s) ds dW u 

Jo Jo Ju Jo Ju 

ft fT ft fT ft 

Y(0,T)- / / a(u,s)dsdu- / / cr(u, s) ds dW u + / f(s,s)ds 

Jo Ju Jo Ju Jo 



Hence 



ft ft fT ft fi 

Y(t,T) = Y(0,T)+ r(s)ds- / a(u,s)dsdu- / a(u, s) ds dW u 

Jo Jo Ju Jo Ju 



t ,-T 



so that 



dY(t,T)= r(t)- [ a(t,s)ds dt - [ a(t,s) 

Jt J Ut 

= [r(t) + A(t, T)\ dt + S(t, T) dW t 



dW t 



and thus 



implies 



as required. 



dp = d{e Y ) = e Y [dY+-d[Y}} 



dP(t, T) 
p(t,T) 



r(t) + A(t,T) + -\\S(t,T)\\ 2 



dt + S(t, T) dW t 



□ 



Example 2.10 Synthetic Money Market Account 

In a bond market, subject to the conditions enumerated before, it is possible to synthetically 
create a locally risk-free bank account. This is accomplished by rolling over just maturing 
bonds. 

Consider a portfolio V which, at any time, consists solely of bonds maturing at time t + dt. 
Suppose that there are nt such bonds in the portfolio, so that 



V t = n t p(t, t + dt) 
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By the self-financing condition, 
dV t = n t dp(t, t + dt) 



= n t p(t, t + dt) 



r(t)+A(t,t + dt) + -\\S(t,t + dt)\y 



Now as dt -► 0, also A(t, t+dt) = - f t +dt a(t, s) ds -► 0, and S(t, t+dt) = - p t +at a{t, s) ds -► 
0. Thus in the limit, 

= r(t)V t dt 

which are just the dynamics of the MMA. 

(Note, however, that the above argument is heuristic in nature: It requires, in any time 
interval, however short, the use of infinitely many types of securities.) 



dt + S(t, t + dt) dW t 



rt+dt 



□ 



In the riskneutral world, discounted bond price processes are martingales, and thus 

p(t,T) =E Q [e-Ji Tr W ds p{T,T)\F t 



e 



St r 0) ds 



In a Brownian world, any equivalent measure is obtained from the objective measure by a 
Girsanov transformation — a consequence of the Martingale Representation Theorem. If 



dF 



e tfudW t -ltf\\u\\ 2 dt 



and if L{t) = Ep[^|jF t ] is the associated likelihood process, then dL t = UfL t dWt- Now if 
p(t,T) = >^p-, then (under Q) 

dp(t,T)=p(t,T)v(t,T) dW t 
(where Wt is a Q-Brownian motion), so that 

dp(t, T) = r(t)p(t, T) dt + p{t, T)v(t, T) dW t 
Hence under P, we have dynamics 

= [r(t) - u(t)v(t, T)] dt + v(t, T) dW t 

i.e. in a Brownian world bond price dynamics are necessarily of the form dp = pm dt+pv dWt- 



3 Modelling the Short Rate 

Short rate models are bond market models where the only explanatory variable is the short 
rate r. This was the earliest approach to bond market models, dating back to the paper by 
vasicek (1977), but short rate models have limited power. Nevertheless, principal component 
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analysis shows that typically 80 - 90% of price variation in the bond market can be explained 
by a single factor, so these models are not wholly devoid of realism. 

When we specify only the short rate, the only exogenously given asset is the MMA At. 
Zero coupon bonds will be regarded not as primitive securities, but as derivatives of the short 
rate. 

Question: Are bond prices uniquely determined by the P-dynamics of the short rate? 

We assume that we live in a Brownian world governed by an objective probability measure 
P, with change driven by a (multidimensional) Brownian motion Wf. We further assume short 
rate dynamics of the form 

dr(t) = (x(t, r) dt + a(t, r) dW t 

i.e. the short rate is an Ito diffusion. 

The answer to the above question is No! 

• The above bond market is clearly incomplete: 

— We are able to execute trading strategies which consist of putting all our money 
in the bank account only. This clearly doesn't give us enough freedom to replicate 
all possible .T^-measurable ransom variables. 

— There is at least one source of randomness, but there are no risky assets. 

— Under any measure , the discounted MMA ^ is a martingale, hence any measure 
equivalent to P, including P itself, is an EMM. The EMM is not unique. 

• If Q ~ P is any equivalent measure, then Q generates an arbitrage-free bond market 
with prices 



p(t,T) =E< 



• In the Black-Scholes model, we also had one source of uncertainty, but there option 
prices are determined by the dynamics of an underlying which is traded. The crucial 
difference here is that the underlying is the short rate, which is not a traded security. 

Nevertheless, bonds of different maturities must satisfy certain internal consistency con- 
ditions in order to exclude arbitrage. For example, if T\ < T2, then p(t,Ti) > p(t,T2), or else 
there will be arbitrage (assuming positive rates). 

If we have d sources of noise (i.e. Wt is a ci-dimensional Brownian motion), then we may 
pick d maturities, and regard the bonds of those maturities as "primitive" securities; bonds of 
all other maturities will be "derivative" . Our market now has as many risky primitive assets 
as securities, and is therefore complete. 



3.1 The Term Structure PDE 

Assume that we have an arbitrage-free bond market, with P-short rate dynamics given by 

dr{t) = n(t, r) dt + a(t, r) dW t 

where Wt is a one- dimensional P-Brownian motion. We restrict to one dimension purely for 
ease of exposition - similar results hold in the multidimensional case. 
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Also assume that the price of a T-bond at time t is given by a sufficiently smooth and 
regular function F: 

p(t,T)=F(t,r(t);T)=F T (t,r) 

By taking two bonds of different maturities S and T, we are able to create a locally riskless 
portfolio. Arbitrage considerations then dictate that the drift of this portfolio is equal to the 
short rate. As usual, this yields a PDE, as we now show. 
First note that by Ito's formula 



dF 1 



if + nFj + l -a 2 Fl 



dt + o-Fj dW t 



(where subscripts denote partial derivatives), so that 



dF T _ 

where 

a T (t,r) = 



a T dt + a T dW t 



F t T + fiF^F + \a 2 F rr 



T aFj 
a (t,r) - 



Consider now a portfolio V consisting of S- and T-bonds with relative weights w s , w T re- 
spectively. Then 



dV s dF s T dF T 



V F s F T 

To eliminate risk, set w a + w T o~ T = 0. Since weights add up to 1, we therefore obtain 



T „S 



Q 0~ rp 

W = ?T W 



a T — a s a T — a s 

Then 

dV _ a s a T - a T a s 
-y ~ a T -a s 

i.e. 

a s a T - a T a s 



a T — a s 



i.e. 

a s (t,r) — r a T (t,r)—r 

Now a T is just the drift of the bond price p(t,T) = F T (t,r), and a T is its volatility. Thus 

a T (t, r) — r 



Market Price of Risk = A 



i.e. all bonds have the same market price of risk A = X(t,r). A is independent of maturity 
(though it may vary over time). 
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Proposition 3.1 (Term Structure PDE) 

In an arbitrage-free one-factor short rate model dr = fi dt + a dWt there is a process X(t, r) 
such that 

a T (t, r) — r 



Market Price of Risk = A 



Hence all bonds satisfy the following PDE 



F T (T, r) = 1 



Proof: Since 



pT 



we have 



* T (t,r) 



Q-F+ 



l a 2 F T - rF 1 

2 rr 11 



cjFJ 



which can easily be manipulated to yield the term structue PDE. 



Using the Feynman-Kac formula, we see that the bond prices are given by 
F T (t,r)=F(t,r;T)=El? x 



p(t,T) = E% 



e J * 



St r 0) ds 



e J t 
P(T,T) 



ft r ( s ) ds 



i.e. 



where 



dr = (fi — a A) ds + a dW s s > t 
r{t) = r 

are the dynamics of r under Q\ Note that, since rt is a Markov process, we have 



E 



t.r 



e -f t T r(s) ds p ( T ^ T ) 



From the fact that 



p(t,T) = E QX \e-f? r ^ ds p(T,T)\F t 



□ 



it follows that each Q x is a risk-neutral measure (i.e. an EMM for the MMA). 

We can also get the risk-neutral short rate dynamics from Girsanov's Theorem: a Girsanov 
transformation which effects the change of measure from real-world to risk-neutral has a 
Girsanov kernel equal to the negative of the market price of risk. Thus — aX is added to 
the drift when we change the measure. Each market price of risk process A gives a different 
risk-neutral measure Q A . 

To summarize: 

• In an arbitrage-free short rate model, all bonds have the same market price of risk, 
regardless of maturity. 
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• Different market prices of risk yield different risk-neutral measures — The bond market 
is not complete. 

• The agents in the market will (implicitly) determine A and thus Q x . 
3.2 Martingale Models of the Short Rate 

We model the short rate directly under a fixed riskneutral measure Q. This is the EMM 
chosen by market participants, and should, in principle, be hidden in the term structure of 
bond prices. By calibrating a short rate model to bond prices, the market price of risk, and 
thus the market EMM, can be determined. This procedure is known as inverting the yield 
curve, and works as follows: 

(1) Choose a short rate model (Ho-Lee, Vasicek, Cox-Ingersoll-Ross, Black-Derman-Toy) 
involving one or more parameters a = (ai, . . . ,a ra ). The Q-dynamics of the short rate 
are given by 

dr(t) = p(t,r(t);a) dt + a(t,r(t); a)dW t 

(2) Solve the term structure PDE. In the risk-neutral world, the market price of risk is A = 0, 
and thus the PDE is 

if + ^Fj + l -o 2 Fl-rF T = Q 
F T (T,r) = 1 

for all maturities T. This yields theoretical bond prices 

p(t,T;a) =F T {t,r t ;a) 

(3) Go to the market, and "observe" the empirical term structure of bond prices {p*(0,T) : 
T > 0}. 

(4) Choose a so that the theoretical prices p(0,T;a) fit the empirical prices p*(0,T) "as 
closely as possible" (where "close" must be defined somehow. For example, one method 
would be to pick maturities T±, . . . , T n and to pick a±, so that 

n 

^(p(0,T;a)-p*(0,T)) 2 
fc=i 

is minimized.) Let a* be this "best" parameter. 

(5) We now have dynamics 

dr{t) = fi(t,r(t);a*) dt + a(t,r(t); a*) dW t 

under the risk-neutral measure. We can also, in principle, observe the real-world dynam- 
ics 

dr(t) = Jldt + a dW t 
Since [i = p, — aX, we now know the market price of risk A, and thus Q = Q A . 
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(6) Ideally, we would like to have 

p{0,T;a*) =p*(0,T) for all T 

However, these are infinitely equations (one for each T), in only finitely many unknowns 
(the ai,..., a n ). This system is over-determined, and the model can not be made to fit 
the initial term structure of bond prices. 

(7) However, if we choose a to be an infinite dimensional vector, rather than a finite dimen- 
sional one, there may be sufficient room to fit the term structure exactly. For example, 
the Ho-Lee model is given by 

dr{t) = 6{t) dt + a dW t 

where a is a constant, and Wt a one-dimensional Brownian motion, here a = {9{t) : t > 0) 
is an infinite -dimensional vector. The Ho-Lee model can be fitted to the empirically 
observed term structure, but this is not obvious a priori. 

(8) Once we've parametrized our model,we can fit other interest rate derivatives. 
3.3 Common Short Rate Models 

The following are common short rate models with just one source of noise: 

• Vasicek: 

dr = {b- ar) dt + a dW t 

where a, b, a are constants. 

• Cox-Ingersoll-Ross: 

dr = (a — br) dt + o\pr dWt 

where a, 6, a are constants. 

• Dothan or Rendlemann-Barter: 

dr = ar dt + ar dWt 

where a, a are constants. 

• Merton: 

dr = a dt + a dWt 

where a, a are constants. 

• Ho-Lee: 

dr = 9{t)dt + a dW t 

where a are constants. 

• Hull- White (extended Vasicek): 

dr = (b(t) - a{t)r) dt + a(t) dW t 
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• Hull-White (extended CIR): 

dr = (b(t) - a(t)r) dt + cr(t)y/r dW t 

• Black-Derman-Toy: 



dr = a(t)r dt + a(t)r dW t 



• Black-Karasinski: 



dr = (a(t)r + b{t)r\nr) dt + a(t)r dW t 
All of the above can be written as 

dr = (ai(t) + a 2 (t)r + a 3 (t)r In r) dt + + /%(t)r)" dWt 

3.4 Term Structure Derivatives 

Consider the general short rate model dr(t) = fi(t, r) dt+cr(t, r) dWt- Suppose that an interest 
rate derivative has a terminal payoff 3>(T, tt) and a dividend rate q(t,rt) over the interval 
[0, T]. The time-t price of the derivative is obtained via an arbitrage argument: Start with a 
portfolio V consisting of one derivative F and — n T-bonds p. Because of the dividends, we 
obtain 

dV = dF -ndp + qdt 
But choosing n = f^/gp will make the portfolio locally riskless, and we obtain 

F t + \a 2 F„ -rF + q _ P t + \a 2 P„. - rP 

dF — dP 

dr dr 

Now the term structure PDE states that 

1 



and thus 



P t + ^a 2 P rr -rP = -(/x - a\)P r 



F t + (n- a\)F r + X -a 2 F„ - rF + q = 



F(t,r T ) = ${T,r T ) 

This is the generalized term structure equation for an interest rate derivative F (where A = 
if we model the short rate in the risk-neutral world). 
The value of the interest rate derivative is clearly 



F(t,rt;T) = E { 



e-/t T ''« du $(r,rr) + J e-J>«%(s,r s ) 



ds 



A trivial generalization of the Feynman-Kac argument shows that the solution of a PDE of the form 



Ft + nF x + ^a 2 F xx -rF + h = F{T, x) = $(T, x) 



is given by 



where 



F(t, x) = E r ' x y e~ f ? rs d3 h(u, X u ) du + e~ $ r ° ds $(T, X T ) 



dX s — n ds + a dW s for t < s < T and X t = x 
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Example 3.2 (a) A call with strike K and expiry r on a discount bond p(t, T) (with T > r) 
has q = and 3>(r, r) = (p(r, T) — To calculate the option price, we first have to 

solve the term structure PDE to get the bond prices, and then once more to price the 
option. 

(b) An interest swap (pay-fixed) can be idealized as a contract paying a divided rate h(t, rt) = 
r t — r*, where r* is the agreed-upon fixed rate (the swap rate at inception). Here $>(t, tt) = 
0, and so 



F(t,r t ) =E C 



J e -K rudu (r,-r*) ds 



Now a floating rate note paying a continuous rate rt must be priced at par = 1 in order 
to avoid arbitrage (Why?), and thus 



e - St ^ du + ^ e - // »•« du^ ^ 



It follows that 



■//r„du ds 



1 -Eq 



l-p(*,T) 



Hence the value of this idealized swap is 

F(t,r t ) = l-p(t,T) -r* J p(t,s)ds 
The swap rate at time t for maturity T sets the value of the swap to zero, and is 

f t p(t, s) ds 

(c) A cap can be idealized as a derivative with zero terminal payoff and a dividend rate 
q(t, rt) = (r t — f) + , where f is the cap rate. 



□ 



We will spend quite a bit of effort pricing options on discount bonds in the next few pages. 
But what about coupon -bearing bonds, which are, after all, more commonly traded in the 
market? Jamshidian's Trick sometimes holds the answer: In a short rate model, a call on a 
coupon bearing bond can be priced as a portfolio of calls on zero coupon bonds, provided 
that the value p(t,T) = p(t,rt;T) of the zero coupon bonds is a strictly decreasing function 
of the short rate. 

Theorem 3.3 Let C K ' T (t,r t ) be the time-t value of a call on a coupon bond B, where r is 
the expiry of the call, and K the strike. Suppose that the coupon bond pays a coupon Y-i at 
date Ti, where r < T± < ■ ■ ■ < Tjy. Let 

K l = p(r,r* ,Ti) 



where r* solves 



B(t, r*) = K 
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Recall that B(r,r) — '}2 i Y i p(T,r,T i ). Since each p(t,r, Tj) is a decreasing function of r, so is B(t,r), which 
implies that r* is unique. Solve numerically for r* , e.g. via bisection method. Then 

C K ' T (t,r) = ^YiC K » T ' T >(t,r) 

i 

where C K,T ' S is the time-t value of a strike K, expiry T call on a zero coupon bond p(t, S) 
(with S > T). 

Proof: The payoff of the call on B is 



(B(t, r T ) - K) + = Y*>(t, ty, Ti) - xj 



Since each p(t,r,T) is decreasing in r, so is B(r,r T ). Let r* be the unique value of r T for 
which the call C expires at the money, i.e. for which 

B{r,r*) = K 

Now define K, L = p(r, r*, Ti). Then 



Now consider two cases: 
Case 1: If r T < r* , then 



J2 Y iP (r, r T , Ti) > £ Y t p(r, r*,Tj = K 

i i 

and 

p(T,r T ,Ti) >p(r,r*,Ti) = Ki 
Thus if C K,T expires in the money, then so does each C Ki ' T,Ti , and 

(^2 Y iP(r, r T , Ti) -Kj =J2 YiP(r, r T , Ti) - K 

= Y J Y i (p{T,r T ,T i )-K i ) 

i 

= Y,^(p(r,r T ,Ti)-Ki) + 

i 

Case 2: If r T > r* , then Y,i Y iP(T,r T ,T$ < K and p(r,r T ,Ti) < K t . Thus if C K < T expires 
out of the money, then so does each C Ki ' T ' Ti , so 



0= ^2Yip(T,r T ,Ti) - Kj 
= J2 Y iiP(r,rr,T i )-K i f 



36 



Affine Term Structure Models 



Hence in either case 



C K ' T (T,r T ) = Y,^C K ^(r,r T ) 



Thus, by the law of one price, 



for all t < t as well. 



3.5 Lognormal Models 

The Dothan, Rendleman-Barter, Black-Derman-Toy and Black-Karasinski all yield lognor- 
mal short rate dynamics. All suffer from the following problem: Let At denote the money 
market account, with dA t = r t A t dt, Aq = 1. Then 

E [e^o rs H wE y (ro+ri) 

for sufficiently small t. Now, since rj is lognormal, define Yt = lnr^. Then we have an 
expectation of the form 



E 



for some normally distributed Z. Now 

/oo 
-ex 



t e Y t 



e e e 



E 



dz 



OO 



as e eZ » e~ z2 / 2 for reasonable values of z. Hence E[ J 4(t)] = oo even if t is small, i.e. the 
bank account, on average, explodes. 
Indeed, it can be shown that 



E 



1 



oo for all i > 



One consequence of this lognormal explosion is that one cannot price Eurodollar futures. 



4 Affine Term Structure Models 
4.1 Mechanics of ATS models 

Definition 4.1 A short rate model is said to possess affine term structure (ATS) if bond 
prices are given by 

p{t, T) = F T (t, r(t)) = e ^)-B{t,T)r{t) 
where A(t,T), B(t,T) are (sufficiently regular) deterministic functions. 



□ 
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Note that not all short rate models are affine term structure models. However, the class of 
affine term structure models is quite well understood: They are those for which both the drift 
and the volatility-squared are affine functions of the short rate. 

For consider a short rate model with risk-neutral dynamics dr(t) = fj,(t, r) dt + a(t, r) dWt 
and suppose that bond prices are of the form p(t,T) = F T (t,r(t)) = e A ( t ' T ^~ B ( t ' T ^ r ( t \ Sub- 
stituting this expression into the term structure PDE 



ft 



T 



1 



t + + -a 2 Fj r -rF T 



F T (T, r) 




1 



we obtain 



A t - l iB + -a 2 B 2 - (1 + B t )r 







Moreover, since p(T,T) = 1, we must have A(T,T) = = B(T,T). 

If we assume that the drift and volatility of the short rate can be expressed in the form 

H(t,r) = a(t)r + P(t) 
a 2 (t,r) = j(t)r + S(t) 



then we obtain 



A t -0B+ -5B 2 



l + B t + aB- \- f B 2 



The lefthand side is independent of r, whereas the righthand side contains r. This can happen 
only if both sides are identically zero, so that we obtain a coupled system of differential 
equations: 

( A t (t,T) = 0(t)B(t,T) - l -5{t)B 2 {t,T) 
A(T,T) = 

B t (t,T) = -a(t)B(t,T) + ^(t)B 2 (t,T) - 1 
[B(T,T)=0 

Note that the bottom equation (a Riccatti equation) does not contain A, and can therefore 
be solved (in principle, although this may be quite hard). The solution can then be plugged 
into the top equation to solve for A. To solve this equation, simply integrate both sides (from 
t to T). 

Thus a short rate model has affine term structure whenever fj,, a are of the form fj,(t, r) = 
a(t)r + (3(t) and o~ 2 (t,r) = 7(£)r + 6(t). The Ho-Lee, Cox-Ingersoll-Ross, Merton, Vasicek 
and Hull- White models all have ATS. The Dothan and Black-Derman-Toy models do not. 



Example 4.2 The Vasicek Model 

Here we have dr t = (b — ar) dt + a dWt, where a,b,a are constants. Thus we have a 
—a, [3 = b, 7 = 0, 5 = a 2 , all constant. 

The system of differential equations that must be solved is therefore 



A t 

[A(T,T) = 



bB - -a 2 B 2 
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j B t = aB-l 
\B(T,T) = 

The bottom equation is a first order linear equation. This can easily be solved: Use e~ at as 
an integrating factor to obtain 



B(t,T) = e 



at 



1 



e -ot + C ( T ) 



and then use B{T,T) = to get C(T) = -\e~ aT ■ Hence 

1 



B(t,T) 

Plug this into the equation for A to obtain 



1 - e -°(r-0 



= b- 
t 

A(T,T) = 
Integrate both sides: 



A t = bB(t,T)-^a 2 B 2 (t,T) 



1 - e~ a ( T -*) 



^ 1 
"2 a ^ 



1 - e~ a{T ~^ 



A(t,T) = A(T,T) - J T A t (s,T) 

f T 1 



ds 



1 - e -< T -°) 


2 

ds — I 


1 - 













2" 7 t a 2 
a 2 B 2 {B _ {T _ t)){ab - iq*) 

-,2 



ds 



4a 

Now that A, i? have been found, bond prices are given by the equation p(t , T) = e A (*' T )~- B ( i ' T ) r M . 

□ 



In order to invert the yield curve in the above example, the parameters a, ba must now be 
chosen so that the model fits empirical (observed) term structure of bond prices {p*(0,T) : 
T > 0} as "closely" as possible. Clearly, however, we have infinitely many bond prices, but 
only three parameters, i.e. the system is highly over-determined, and therefore we cannot 
generally choose a,b,a such that e A (o,T)-B(o,T)r _ p*(o,T), i.e. the model cannot be made 
to fit the observed term structure exactly (unless we are astoundingly fortunate). The Vasicek 
model is able to fit, exactly, just 3 bonds. 

Example 4.3 Cox-Ingersoll-Ross model The risk-neutral short rate dynamics assumed are 

dr t = (b — ar t ) dt + a^pF t dWt, a, b, a, ro > 

This is mean reverting (to b/a). Since the volatility term a^/rl tends to zero as rt — > 
(which is consistent with observation), positive rates are assured (which is also consistent with 
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observation) . Postulating p(t, T) = e A ( t ' T ^ B ( t ' T ) rt , we quickly determine, by substituting into 
the term structure PDE, that 

B t = aB + ^or 2 B 2 - 1 B(T, T) = 
A t = bB A{T, T) = 

To solve the Riccati equation for B, we try a solution of the form 

X{t) 



B(t,T) 



Then 



Bt = 



cX(t) + d 
X t cXX t 



cX + d (cX + d) 2 
and hence, substituting into the equation for Bt, we see that 

-dX t + X 2 {ac + ^a 2 - c 2 ) + X{ad - 2cd) - d 2 = X{T) = 

Choose c to ensure that a + \a 2 — c 2 = 0, i.e. c = \{a + \/a 2 + 2cr 2 ). We then have a order 
linear differential equation 

X t + kI = — d where k = — a + 2c = \J a 2 + 2a 2 
Since X(T) = 0,we see that 



X(t) = -[e K(T_t) -l] 



Hence 



B(t,r) 



cX(t) + d 

e K(T-t) _ 1 

i(K + a)(e K ( T -*) - 1) + k 



2k + (a + K)(e K ( T -*) - 1) 
Then T) is obtained by integrating: 

cT 



2( e *(T-t) _ i) , 

— — t^—t where k = v a 2 + 2a 2 



The solution is 



A(t,T) = A(T,T) - f A t (s,T) ds = -b f B(s,T) 
Jt Jt 

2 Ke K a + K )( T -*) 



ds 



2b 

A(t, T) = -2 la 



2K + (a + K)(e K ( T "*) - 1) 



as can be verified by differentiation. 



□ 
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Example 4.4 Ho-Lee Model 

We are given risk-neutral short rate dynamics dr{t) = 6{t) dt + a dWt, where 9{t) is deter- 
ministic and a a constant. The model has ATS with a = 0, = 9, 7 = 0, 5 = a 2 . This leads 
to two differential equations. The first is 

f B t = -1 
[ B(T, T) = 

which has solution B(t, T) = T — t (as can be seen by integrating both sides from t to T). 
The second DE is 

A t = e(t)B(t,T)-^a 2 B\t,T) 

= e{t){T-t)-\a\T-tf 

[A(T,T)=0 
Integrating both sides from t to T yields 

cT 



A(t,T) = - [ 6{s){T-s)ds+ l -<T 2 {T-tf 
Jt o 



We now choose the function 9(t) so as to fit the initial term structure of bond prices {p*(0, T) : 
T > 0}, or, equivalently, the observed term structure of (instantaneous) forward rates 
{/* (0,T) : T > 0}. 

Recall that f*(0,T) = - 81n ff' T) . With affine term structure, we have p*(0,T) = 

e A(0,T)-B(0,T)r 0) g0 



lnp*(0 

Differentiating with respect to T, we see that 

pT 



, T) = - f 9{s){T -s)ds + \o 2 T* - r T 
Jo o 



r(0,T) = J 6{s)ds- l -a 2 T 2 + r Q 



Differentiating once more with respect to T, we obtain 



df*(0,T) 
dT 



6{T) - a 2 T 



and thus we have found 8: 



8(t) = f*(0,t)+a 2 t 



We can use this to calculate A(t,T): 

fT 1 

A(t,T) = J (f*(0,s)+a 2 s)( S -T)ds + -a 2 (T-t) 
= f*(0,s)(s-T)\J- ^ f*(0,s) ds + a 
= f*(0,t)(T-t) + J 



s 3 s 2 Tl T 



a 2 (T-ty 



r(0,t)(T-t) + ln 



dlnp*(0,s) , 1 o /m x2 

P*(0,T) \ 1 
p*(0,t) 



a 2 t(T-tf 
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Using the fact that the Ho- Lee model has ATS, we see that bond prices are given by 

T) = exp ( r (0 ' t){T ~ t] ~ \ aH{J ~ 1)2 ~ {T ~ t)r{t) ) 

where ^ 

r(t) = r + [ 0(s) ds+ [ a dW s 
Jo Jo 

= r + /* (0, t) - f (0, 0) + l -a 2 t 2 + aW t 
= r^,t) + l -a 2 t 2 + aW t 

because ro = /*(0, 0). It follows that E[rt] — > oo (under the riskneutral measure). This is 
clearly a flaw in the model. 

Since the short rate is Gaussian, future bond prices are lognormally distributed under the 
risk-neutral measure. In particular, there is a non-zero probability that a bond will, at some 
future date, trade above par (i.e. that interest rates become negative). This is clearly another 
flaw in the model. 

Now that we've calculated the evolution of future bond prices and rates, let's have a look 
at future forward rates. Since f(t,T) = _ 81n P(f> T ) ^ we gee ^ na ^. 

f{t, T) = f* (0, T) - r (0, t) + a 2 t(T -t)+r t 
= f*(0,T) + a 2 t(T-h) + aW t 

using the expression for r t obtained earlier. Note that f(t, t) = r t . 

Now if we fix t > 0, we see that E[/(£, T)] -^ooasT^oo. Indeed, for large values of T, 
f(t, T) ~ kT. Thus even if the initial forward curve is bounded above, it will be unbounded 
an instant later. This is another flaw in the Ho-Lee model. 

□ 

Example 4.5 The Hull- White (extended Vasicek) Model 
Consider the short rate model with risk-neutral dynamics 

dr t = {b(t) - ar t ) dt + a dW t 

where b{t) is deterministic, a, a are constants and W% is a one-dimensional Brownian motion. 
This is clearly an affine term structure model drt = (a(t)rt + (3(t)) dt + y/j(t)rt + 5(t) dWt, 
with a(t) = -a,P(t) = b,j(t) = and S(t) = a 2 . Substituting p(t,T) = e A(t,T)-B(t,T)r t into 
the term structure PDE yields 

B t (t,T) = aB(t,T) - 1 
B(T,T) = 

and 

A t (t,T) = b(t)B(t,T) -^a 2 B 2 (t,T) 
A(T, T) = 
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Hence 

B(t, T) = i(l - e - a ( T "*)) 

A(t, T) = [ -b(u)B(u, T) + \a 2 B 2 {u, T) du 
Jt 2 

Fitting the initial term structure of bond prices is equivalent to fitting the initial term struc- 
ture of forward rates. The latter is more convenient. Now since /(0,T) = — dllL g^ l ' T ^ = 
-A T (0,T) + B T (0,T)r , and since B T (t,T) = e - a ( T -*), we observe 

/(0,T)= / b{u)B T {u,T) + a 2 B T {u,T)B(u,T) du + B T (0,T)r 
Jo 

= [ T b(u)e- a ^ du - ^(1 - e-< T ~^f + e- aT r 
Jo 2a 

We side-step the computation of these integrals using the following trick: Define 

x(t) = e- at r + f b(u)e- a ^ du 
Jo 



2 

a '- -at\2 



Note that 



X '(T) = -ar e~ aT + b{T) -a I b{u)e~ a ^' u ^ du 

Jo 

= -ax{T) + b(T) 
Now /(0,T) = x(T) - y(T), and so 

6(T) = x + ax 

= / T (0,T)+y / (T)+ax(T) 

= / T (0,T)+y / (T)=a[/(0,T)+y(T)] 

Thus, noting that y(t) = |J(l-e~ at ) 2 = \o 2 B 2 (<d,t) and thus that y'(t) = ^(l- e - at )e- at = 
a 2 B(0, t)B T (0, t), we obtain 

b(t) = f*(0, t) + a 2 B(0, t)B T (0, t) + a[/*(0, t) + ^ 2 # 2 (0, t)} 

= / * ( 0,t) +a r(0,t) + ^[l-e- 2 i 

This is the function b{t) which will fit forward rates to the observed term structure {/*(0, T) : 
T > 0}. 

Since we now know b(t) we can calculate A(t,T): 

A(t,T)= [ -b(u)B(u,T) + \o 2 B 2 (u,T) du 
Jt 2 

Now note that b{u) = x'{u) + ax(u) = e~ au d X J^ U so that 
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b(u)B(u,T) du 



,dxe au 
du 



(1 - e - a{T ~ u) ) du 



1 f T 

= - / (e' au - e~ aT ) d{xe° 
a Jt 



[x{u)e au {e 



mi i , - n,u g— aT^j 



T _ 1 

* a 



xe au ■ -ae au du 



-x(t)(l - e~< T -V) + f x(u) du 
1 Jt 

B(t,T) + J - 



f(0,t) + Y B 2 (0,t) 



dpjO, u) a 2 2 



-/(O, t)B(t, T) - In - ^5 2 (0, t)B(t, T) 



T a 2 



B 2 (0,u) du 



J B 2 (u, T) - B 2 (0, u) du + B 2 (0, t)B(t, T) 



Hence 

A(t,T) = f(0,t)B(t,T) + ln^^ + ^ 

Now, after a few lines of manipulation, 

f B 2 {u,T) - B 2 {0,u) du + B 2 {0,t)B{t,T) = -—B 2 {t,T)(l 
Jt 2a 

as you can easily check, substituting B(t,T) = \{l - e" a ( T -*)). Thus 

A(t, T) = /(0, t)B(t, T) + In - ^B 2 (t, T) (1 - e - M ) 

Substituting p(t, T) = e A{t,T)-B{t,T)n we obtain: 



^-2at\ 



p(t,T) 



P(0» T ) o /(0,i)fi(( ) T)^B 2 ((,T)(i-e- J « l )-B(i,T)r l 



p(0,t) 



□ 



We have thus found the following bond prices: 

Theorem 4.6 (a) In the Ho-Lee model, bond prices (fitted to the initial term structure) are 
given by 



p(t,T) = 



P(0,T) 
p(0,t) 



exp /(0, t)(T - t) - -a 2 t(T - t) 2 - (T - t)r(t) 



(b) In the Hull-White (extended Vasicek) model, bond prices (fitted to the initial term struc- 
ture) are given by 



p(t,T) 



p(0,T) 



p(0,t 

where B(t, T) = \{l - e^"*) 



• ./;()./)/>•;/. /•) °^B 2 - < 2 "')- B(i.T)n 



□ 
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4.2 Bond Options 

In the chapter on changes of numeraire, we obtained the following general option formula: 
The price of a call C with strike K and maturity T on an underlying S is given by 

C = S Qs(S T > K) - Kp(0, T)Q T (S T > K) 

where Qs,Q are the EMM's associated with numeraires St,p(t,T) respectively. 

In order to use this formula, and to get Black-Scholes type solutions to option pricing 
problems, we assumed that the volatility of the securities is deterministic, and then obtained 

Theorem 4.7 // St = p (fy) *s an Ho process of the form = fi(t) dt + a(t) ■ dWt, and if 
a(t) is deterministic, then the value of a call C with maturity strike K and T on underlying 
security S is given by 

C = S N(d 1 )-Kp(0,T)N(d 2 ) 

where 

0~av 

di 

Put-call parity yields 
for the price of a corresponding put. 

□ 



m Kp(0,T) + 2 a av 1 r- 

= ^ — ^7= d 2 = di - o-avVT 

P = -SoNf-d!) + Kp(0, T)N(-d 2 ) 



We can now use this theorem to price bond options. 



Example 4.8 Bond Options in the Ho-Lee Model 

Consider a European call option C with strike K and maturity T on a discount bond p(t, S) 
(where S > T). In the Ho-Lee model, with risk-neutral dynamics dr t = 6(t) dt + a dWt, 
bond prices have dynamics 

= rdt-a{T- t) dW t 
p(t,T) V ; 

The drift term is r, because bond prices have drift r under the risk-neutral measure, just 
like all other traded securities. The volatility is obtained from the affine term structure: 
p(t,T) = e Mt,T)-B(t,T) rt > an j we f oun j t h a t B(t,T) = T -t (and we don't care about the 
value of A(t, T) right now.) Thus the bond volatilities are deterministic: p(t, T) has volatility 
—a(T — t) and p(t, S) has volatility —a(S — t). Now the underlying security is p(t, S), and 
p(t, S) = ^'t] has deterministic (indeed, constant) volatility —a(S—t)+a(T—t) = —a(S—T). 
This is because the volatility of a ratio of two assets is just the difference of their volatilities. 

It follows that p(t, S) is lognormally distributed, and that lnp(t, S) has variance c^ v T = 
J T <j 2 (T - S) 2 dt = a 2 {S - T) 2 T. It follows that the price of the call is 



C = p(0, S)N(d 1 ) - Kp(0, T)N(d 2 ) 
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where 

ui = 

(j(S-T)VT 

d 2 = di -a{S-T)yff 

□ 

Example 4.9 Bond Options in the Hull- White (extended Vasicek) Model 
We tackle once more the problem of pricing a call with strike K and maturity T on a zero 
coupon bond p(0, S), where S > T. It ought to be clear from the analysis of bond options in 
the Ho-Lee model that we need mainly to find the volatility of the bonds p(t,T). Now, as 
for the Ho-Lee model, the riskneutral dynamics of p(t, T) are 

^2^1 = r dt - B(t,T)a dW t 

so that the volatility of p(t,T) is —-(1 — e" a ^ _ *- ) ). The asset ratio pt = p(t, S) /p(t,T) 
therefore has volatility ^ (e~ aS — e~ aT )e at at time t. Thus the average volatility-squared is 

a\ v T = f T °^{e~ aS - e~ aT ) 2 e 2at dt = ^(1 - e- a ^) 2 (l - e~ 2aT ) 
Jo a 2a 

We now find that the value of the call is simply 

p(0,S)N(d 1 ) - Kp(0,T)N(d 2 ) 

where 

, 111 Kp(0,T) ^ 2 a av 1 r- 
dl = —j= d 2 = dl- (TavVT 

□ 



5 The Heath— Jarrow— Morton Framework 
5.1 The Set-Up 

Up till now, we have studied interest rate models in which the short rate is the only explana- 
tory variable. Such an approach has many obvious advantages: 

• Specifying r as the solution of an SDE allows us to use Markov theory, which leads to 
PDE's (e.g., via the Feynman-Kac theorem, or the Kolmogorov forward and backward 
equations) that can be solved; 

• If we're lucky, we can obtain analytical formulas for bond prices and bond option prices 
(as we did for the Ho-Lee and Hull-White (extended Vasicek) models. 

However, the short rate modelling approach has some obvious disadvantages as well: 

• It is unreasonable to regard the short rate as the only explanatory variable — it is 
difficult to incorporate views about different times in the future; 
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• It can be quite difficult to fit a realistic volatility structure; 

• In order for the model to have even a remote chance of being correct, it is necessary to 
invert the yield curve (i.e. to fit the model to the initial term structure of bond prices). 
This can be quite difficult as well. 

The Heath-Jarrow-Morton (HJM) approach circumvents some of these difficulties by 
specifying dynamics for the entire (uncountable) family of forward rates. For a fixed T > 0, 
assume that the forward rate f(t,T) has "real-world" dynamics 

df(t, T) = a(t, T) dt + a(t, T) dW t T > 0, <t <T 

where Wt is a finite-dimensional Brownian motion under the real world measure P, and 
a(t, T) and a(t, T) are adapted (and sufficiently regular to ensure that most of the operations 
below are permissible. For example, it is often necessary to assume that a(t, T) is jointly 
measurable in the t- and T-variables.) 

Thus we have infinitely many SDE's, one for each maturity T. Each such SDE has an 
initial condition, namely /(0,T) = /*(0,T), where /*(0,T) is the observed term structure, 
the advantage of this approach is that the initial term structure is fitted automatically — it is 
an initial condition! — so that inverting the yield curve becomes unnecessary. It is also easier 
to incorporate views about different maturities, because we have many different SDE's. (The 
disadvantage, of course, is that we have many, many SDE's.) These are still manageable, 
because we assume that the bond market is driven by finitely many sources of noise. But this 
leads to another difficulty: 

Remarks 5.1 Given a(t, T), a(t, T) and {/*(0,T) : T > 0}, we can solve the SDE's for the 
forward rate, so that we have specified the entire term structure {f*(t, T) : T > 0,0 < t < T} 
at all times and all maturities, and thus the entire term structure of bond prices 

p(t, T) = e~ £ f{t ' u) du 

Since we have only finitely many sources of noise, and infinitely many traded assets, there is 
a possibility of arbitrage in the bond market, unless the bond prices are inter-related in a 
specific way (which amounts to all bond prices having the same market price of risk, for all 
source of noise). This will impose conditions on the functions a and a. 

□ 

Remarks 5.2 HJM is not a model, but a framework of models for the bond market; short 
rate models are another such framework. But whereas short rate models are generally ltd 
diffusions, and thus Markov processes, we can easily let a and a depend on past history. HJM 
models therefore need not be Markov models. (Of course, short rate models do not really need 
to be Markov either, but then their dynamics cannot be given by diffusions. We shall discuss 
the relationship between short rate and HJM models in the next section.) 

□ 

For a market model (driven by Brownian motions) with only finitely many securities, we know 
that the model is arbitrage-free if and only if we can construct a risk-neutral measure, and 
complete if that measure is unique. Equivalently, the market is complete if an only if there 
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are as many traded risky securities as Brownian motions, subject to some conditions which 
ensure that the traded securities, are, in some sense , independent (where "independent" 
is meant in the sense of linear algebra, and not probability). The fact that there are only 
finitely many sources of noise, but infinitely many traded assets, means that the market is 
"over-complete", i.e. that there may be many ways of replicating a security. Unless all such 
replicating portfolios have the same price, there will be arbitrage. In practice, all securities 
must have the same market price of risk. If that's the case, we can construct a riskneutral 
measure (via a Girsanov transformation), which implies that the market is arbitrage-free. 
The arbitrage theory we've developed thus far only applies to markets with just finitely 
many traded securities, and it isn't at all clear that the impossibility of arbitrage implies 
the existence of a riskneutral measure (i.e. a measure under which all uncountably many 
zero coupon bond prices, when discounted, become martingales). We can, however, construct 
riskneutral measures for any finite subset of zero coupon bonds. Nevertheless, it is highly 
desirable to have a single riskneutral measure for all bonds simultaneously (because prices 
of securities are then just expected discounted payoffs, where the expectation is taken w.r.t. 
the riskneutral measure) . We will therefore try to impose a strong form of the no-arbitrage 
condition: The existence of a riskneutral measure for all bonds. 

To enable us to construct such a riskneutral measure, there must be relationships between 
a(t,T) and a(t,T) that must hold if the HJM model is to be arbitrage — free: 

Proposition 5.3 Assume that the bond market is arbitrage-free in the strong sense, i.e. 
assume that there is a risk-neutral measure for bonds of all maturities. Then there is a 
(multidimensional) process X(t) such that, for all maturities T, 



If we use a Girsanov transformation with kernel —A to change to a new measure Q, then new 
dynamics of p(t, T) are 




Proof: Recall that 



dp(t,T) 
p(t,T) 



r(t) + A(t, T) + -\\S(t, T) || 2 dt+ S(t, T) dW t 



where Wt is a P-Brownian motion. Here 




dp(t,T) 
p(t,T) 



r{t) + A(t,T) + ^\\S{t,T)\\ 2 - S(t,T)X(t) dt + S(t,T) dW t 



where Wt is a Q-Brownian motion. For Q to be a riskneutral measure, each p(t, T) must have 
drift r(t), i.e. 

A(t,T) + 1 -\\S(t,T)\\'-a(t,T)X(t) = 
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This shows that A is just the market price of risk of p(t, T) at time t, for all T: All bonds 
have the same market price of risk. 

Differentiating this equation with respect to T yields 

-a(t,T) + a(t,T) J a{t, s) tr ds + a(t,T)X(t) = 

□ 

Suppose that we have an HJM model driven by d sources of noise, so that each a(t, T) is a 
(i-dimensional row vector a = (a±, . . . , ad), and A = (Ai, . . . , Xd) tr is a d-dimensional column 
vector. We then have 

d „t d 

a(t,T) = Y,°i(t,T) o-i(t,s)ds + Y,Vi(t,T)\i(t) (*) 
i=i Jt i=i 

If we take a and a as given, we can try and solve for A. We then have uncountably many 
equations in just d unknowns Ai(i), . . . , \d(t) — one equation for each T. Thus a, a cannot 
be specified arbitrarily. What we can do is 

• Specify the volatility surface a(t,T). 

• Choose d benchmark maturities T\, . . . , Td and specify a(t, Ti), . . . , a(t, Td). 

• Solve the system (*) of d equations for the d unknowns Ai(t), . . . , Xd(t)- 
All the other a(t,T) (for T^a bench mark maturity) are now given by (*). 

5.2 Martingale Modelling 

As for short rate models, it is often convenient to bypass the necessity of estimating the 
market price of risk, and to model directly under the risk-neutral measure Q. i.e. we write 

df(t,T) = a(t,T)dt + a(t,T)dW t T > 0, < t < Tf(0, T) =/*(0,T) 

where Wt is a Q-Brownian motion. Under Q, the market price of risk is A = 0, so we obtain: 

Proposition 5.4 (HJM Drift Conditions) 

The riskneutral dynamics of forward rates satisfy the following conditions: 

a (t,T) =a{t,T) J a(t,s) tr ds 

□ 

Thus in the riskneutral world, the drifts a(t, T) are completely determined by the volatility 
surface a(t, T). To create an HJM model, therefore, just follow the following steps: 

• Estimate (or otherwise specify) a volatility surface a(t,T). 

• Calculate the drifts a(t,T) = a(t,T) a(t,s) tr ds. 
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• Observe the term structure of forward rates {/*(0, T) : T > 0}. This involves building 
a yield curve for all maturities. 

• Integrate: 

f(t,T)=f(0,T)+ f a(u,T) du+ f a(u,T) dW u 
Jo Jo 

• Compute bond prices p(t, T) = e~ h ^(*' s ) ds and the prices of other interest rate deriva- 
tives. 



5.3 Examples and Applications 

Example 5.5 We consider here the simplest possible HJM model: We have only one source 
of noise, and put a(t, T) = a = constant for all t, T. By the HJM drift conditions, we see 
that 

f T 

a{t,T)=aJ a ds = a 2 (T - t) 

under the riskneutral measure. Hence the riskneutral dynamics of forward rates are 

df(t, T) = o 2 (T -t)dt + a dW t 
/(0,T) = /*(0,T) 

Integrate this to obtain 

f(t, T) = f* (0, T) + a 2 t(T - * ) + aW t so that 
r(t)=f(0,t) + ^a 2 t 2 +aW t 
and thus the short rate dynamics are given by 

dr t = [/t(0, t) + a 2 t] dt + a dW t 

These short rate dynamics should be familiar: We've obtained the Ho-Lee model fitted to 
the initial term structure! Note that we didn't have to do the actual fitting — in the HJM 
framework, fitting is automatic. 

□ 

Thus the Ho-Lee model is (equivalent to) the simplest HJM model. 

Example 5.6 Can the Hull-White (extended Vasicek) model be recast in the HJM frame- 
work? 

Indeed it can. The Hull- White model drt = {b{t) — art) dt + a dWt is an affine term structure 
model, with bond prices p(t,T) = e A(t,T)-B(t,T)r t _ Hence f(t,T) = -A T (t,T) + B T (t,T)r t , 
which means 

df(t, T) = [■} dt + B T (t, T)a dW t 

where we haven't bothered to calculate the coefficient of the dt-teim (which is, of course, just 
a(t,T)). But for the Hull-White model, it was easy to calculate B(t,T) = ^[1 - e~ a{T -% so 
that B T (t,T) = e~ a( - T - t l It follows that 

df(t, T) = a{t, T) dt + ae~< T -^ dW t 
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Thus a(t,T) = ae a ( T *) We can now use the HJM drift conditions to calculate a(t,T) = 

a(t,T)tfa(t,s) ds = £[ e -a(T-t) _ e -2a(T-t)]_ 

To verify that the above model leads to the Hull-White model, recall that the short rate 
dynamics can be deduced from the forward rate dynamics as follows: 

dr t = [fr(t, t) + a(t, t)\ dt + a(t, t) dW t 

Now a(t, t) = a, and a(t, t) = 0. Finally, 

f(t,T)=f(0,T)+ f a(u,T) du+ f a(u,T) dW u 
Jo Jo 

which implies that 

r(t) = 0(t)+ f ae-^-'^ dW u 
Jo 

for some function &(t), and hence that 

dr{t) = @'(t) dt-(aj ae-< 1 -^ dW^j dt + adW t 

= [e'(t) - a(r(t) - 9(t))] dt + a dW t 
= [b(t) - ar(t)} dt + a dW t 

Moreover, b(t) = @'(t) + a@(t), and @(t) = f(0,t) + J^a(u,t) du = f(0,t) + ^(1 - e~ at ) 2 . 
This is exactly the value of b(t) which we obtained for the Hull-White model fitted to the 
initial term structure. 

□ 

Remarks 5.7 The above example suggests a simple mechanism for turning a fitted affine 
term structure model dr t = fi r dt + a r dWt into an HJM model: 

• If p(t,T) = e A{t,T)-B{t,T)v t j golve the ( Riccatt i) ODE for B{t,T). 

• Then the HJM volatility surface is given by a{t,T) = B T (t,T)a r dW t . 

• The HJM drift conditions now specify a(t, T) as well. 

□ 

Example 5.8 We consider a model with two sources of noise W/, W t 2 and a volatility surface 

a(t,T) = (a 1 ,a 2 e- a ^) 
where a±,a2,a are positive constants. The HJM drift conditions dictate that 



2 

a(t,T) = af(T-t) + ^ 

a 

Integrating the forward rate dynamics, we see 



e -a(T-t) _ e -2a(T-t) 



f(t, T) = /(0, T) + a\t (t - + ^ [2e- aT {l - e at ) - e~ 2aT (l - e™)] 

+ a 1 W^ + a 2 f'e-^-^ dW 2 u 
Jo 
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Thus 



= f(0,t) + + || [2(e-* - 1) - (e- 2at - 1))] 

+ CT!^ 1 + (72 /' e- a (*~«) aW 2 




Thus the short rate is 



a Gaussian process, and 



dr t 




Jo 



[e'(t) - a(r t - 9(t) - ai W/)] + (7i dV^ 1 + a 2 dW t 2 
[&(*) - ar t - aaiW/] dt + ai rfW^ 1 + a 2 dW? 



This is not the form of one of our standard short rate models, because of the explicit presence 
of Wl in the drift. 



6 Market Models: Preliminaries 

The HJM aproach studies the entire term structure of instantaneous forward rates {f(t,T) : 
t < T}, with considerable success, as we have seen. Nevertheless, forward rates for only a 
few maturities are available in the market, so the forward rate curve, like the instantaneous 
short rate, is a purely mathematical entity, a mathematical idealization. Market models, on 
the other hand, model observable (i.e. market-quoted) rates rather than idealized entities, 
and thus simple, discrete rates. 

The London Interbank Offer Rates (LIBOR), for example, are quoted for different ma- 
turities (3-month, 6-month, etc.) and also for different currencies. These LIBOR spot rates 
imply LIBOR forward rates using an arbitrage argument. New LIBOR quotes are available 
daily. Swap rates (the fair rates for interest rate swaps) are another example of discrete 
market-quoted rates. The market model approach to interest rates dates back to Miltersen, 
Sandmann and Sondermann (1997), Brace, Gatarek and Musiela (1997) and Jamshidian 
(also 1997). Several other approaches now exist, due to Hunt and Kennedy, and Musiela 
and Rutkowski, amongst others. It remains one of the most intensively researched areas of 
financial mathematics. 

6.1 Black's Models 

Black's model has long been the industry-standard model used by traders to price a variety 
of European-style options, including interest rate options, such as caps, floors, and swaptions. 
It is essentially a minor variation on the Black-Scholes formula, as we shall shortly see. 
Nevertheless, the suitability and adequacy of Black's model has often been questioned by 
academics, particularly in the arena of interest rate options. 

Consider a European call option C with strike K and maturity T on some market variable 
X. X need not be a traded instrument — it could also be a market-quoted interest rate, 
for example. The main assumption is that Xt is lognormally distributed in the riskneutral 



□ 
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world.. Thus we make no assumptions on the distribution of the process (X t )t in general, but 
just on the value of X at the expiry of the option. We further define the "volatility" of Xt 
to be a non-negative number a satisfying 

variance of In Xt = o 2 T 

Let Q be a riskneutral measure. Then the t = 0-value of the call is 

Co = E Q \e~ n d \X T - K + ) 

Black uses two approximations to determine the value of Cq: 
• Approximate 



e~fo ^ d \X T - K) + \ « P(0,T)E Q [(X T - K)+] 

i.e. discount outside the expectation operator. 

• Now because Xt is lognormal under Q, we know that 

E Q [(X T - K) + ] = E[X T ]iV(di) - KN{d 2 ) where 
In M*r\ , l a 2 T 

d! = - K 2 



Approximate 



aVT 
d 2 = di- aVf 

E,q[Xt] = forward price/rate of X = Fq 



i.e. approximate the expectation by the forward price/rate. 

Since the forward price of X at time T for time T is just itself (i.e. Ft = Xt) , this can be 
interpreted as saying that the forward rate process has zero drift, i.e. is a Q-martingale. 

Thus, using these two approximations, we obtain Black's model for a call on X: 

C = P(0, T)[F N(d 1 ) - KN(d 2 )} where 

In io , i a 2 T 
di = K + 2 

d 2 = di- aVT 

A similar formula is obtained for puts, using put-call parity. 

If payments are based on a variable Xt, but only received at some later date T*, then dis- 
counting must be done from time T* rather than from time T. Black's model then generalizes 
to give call prices 

C = P(0,T*)[F N(d 1 )-KN(d 2 )] where 
lnf + WT 

ii 



A = ~ K 



o-Vt 

d 2 = di- aVf 



Interest Rate Modelling 



53 



where To is still the T-forward value of X at time t = 0. The appropriate generalized Black 
formula for put options follows once again by put-call parity. 

Now it ought to be clear Black's model has several flaws. Firstly, it cannot be appropriate 
to use the first approximation when Xt depends on interest rates, as it amounts to saying 
that 



E Q e'fo ^ d \X T -Ky 



- / r t dt 



E Q [(X T - K)+] 



which is close to asserting that r and Xt are independent. That's a dangerous assumption 
if X happens to be an interest rate derivative! There is no justification for the second 
approximation either. The expected value of Xt under the riskneutral measure is its futures 
price, whereas the forward price is the expected value of Xt under the T-forward riskneutral 
measure. These measures are not the same if interest rates are stochastic. 

In spite of these flaws, Black's model remains heavily used — the industry standard. The 
method can be justified, provided that the relevant variable is taken to be lognormal under 
a different measure, associated with a different numeraire. We shall give several examples of 
this below. Review material on changes of measure and numeraire may be found in the next 
subsection. 



Example 6.1 Bond Options: Lognormal prices 

We consider a call C with strike K and maturity T on a coupon bearing bond B. We assume 
that the bond price at time T is lognormally distributed (under the riskneutral measure) , and 
that In Bt has variance a 2 T. This "volatility" a is obtained from historical data (or implied 
by other market variables). 1 

The T-forward bond price To is simply the fair price which sets the value of a forward 
contract on B equal to zero. A simple arbitrage argument shows 

_ Bq-D 
P(0,T) 

where -Bo is the current value of the bond, T(0, T) is the discount bond maturing at time T, 
and D is the present value of all coupons (dividends) paid out during the life of the option. 
Thus Black's model determines 

C = (B - D)N(di) - KP(0, T)N(d 2 ) where 

i n b q -d I ±rr 2 T 
, 111 KP(0,T) + 2° 1 

d\ = — 



d 2 = di- oVT 



The above call price is an approximation under the assumption that Bt is lognormal under 
Q, but exact if we assume lognormality of Bt under the T-forward riskneutral measure Q T . 
Of course, bond put options can be evaluated by put-call parity. 

□ 

1 In practice, yield volatilities are often obtained. If <j y is the volatility of the yield (i.e. if a y T is the 
standard deviation of the logarithm of the forward yield kiy T ), then (with D* — duration) we have ^ w 

— D* Ay = —D*yo^, i.e. A(lnB) w — D*j/oA(ln y). Thus the variance of InB is approximately (D*j/o) 2 x the 
variance of lny, i.e. ob ~ D*y a y . 
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Example 6.2 Caps: Lognormal LIBOR Rates 

An interest rate cap is an option-like contract which protects the holder against a floating 
interest rate moving too high. Each cap is a portfolio of caplets, each for a certain future 
time interval. A caplet is essentially a call option on the floating rate, given a certain cap 
rate as strike, based on a given notional amount. Consider, for example, a five-year cap, on 
a notional amount A, with cap rate R and semiannual resets based on 6-month LIBOR. This 
is a portfolio of 10 caplets. The reset dates To = 0, T\ = 0.5, T 2 = 1, . . . T\$ = 5 are referred to 
as the tenor structure of the cap. The n th caplet protects the holder against 6-month LIBOR 
rising above R over the period [T ra _i,T„]. It is a call option with strike R on the 6-month 
spot LIBOR L(T n _i) at time T„_i, and will have the following payoff at time T n : 

Payoff of n th cap = A6 n {L(T n -{) - R) + where S n = T n - T n _i 

(This is a payment-in-arrears cap. The first caplet is generally excluded from the cap, because 
there is no uncertainty about the spot LIBOR L(Tq).) 

To price the n th caplet using Black's model, we assume that the future spot LIBOR 
L(T n _i) is lognormally distributed, with volatility a n -\. The t = 0-forward LIBOR rate (i.e. 
the Fq of Black's model) for the period [T„_i,T n ] is given by 

L(0,T n _ 1 )= P(0 ' T "- l) - P(0 ' T " ) 



6 n P(0,T n ) 

(In this notation, the future spot rate, L(T n _i), is just L(T n _i, T„_i).) 
Hence the t — 0-value of the n th caplet is 

C„(0) = A5 n P(0,T n ) [L(0,T n _i)JV(di, n _i) - RN(d 2 , n -i)} 

«l,n-l - /Tp 

dl,n-\ = ^l,n-l — Cn-l\/^n-l 

The price of the cap is therefore the sum of the prices of the caplets (though, as we have 
mentioned, the first cap is often excluded, i.e. C±(0) is set to zero). 

The above price for a cap is an approximation, assuming that each future LIBOR spot rate 
L(T n ) is lognormal under the riskneutral measure Q. The formula for each caplet is exact, 
however, if it is assumed that L(T n ) is lognormal under the T n+ i-forward measure. For then 
indeed 



C„(0) _ F 



P(0,T n ) 



A6 n [L(T n -i) - R} + ] 



P(T n ,T n ) 

which justifies the first approximation used in Black's model (i.e. discounting outside the 
expectation). Moreover, the second approximation is exact, i.e. the forward LIBOR rate 
L(0, T„_i) is exactly equal to the expected value of the spot rate, but under the forward 
riskneutral measure: L(0, T„_i) = Eqt„ [L(T n _i)]. To see this, note that a long forward rate 
agreement F, initiated at time t = for period [T„_i,T n ], will have initial value Fq = 0, and 
terminal value Fx n = <5 n [L(T n _i) — L(0,T n _i)]. Hence 



Ft 

1 n 



P(T n , T n 

which yields the required result (because L(0, T n _i) is a known constant). 
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So in order for the Black price of a cap to be accurate, we must simultaneously assume 
that each L(T n ) is lognormal under Q Tn + 1 . This seems difficult to justify theoretically. One 
of the achievements of LIBOR market models is that they provide a framework under which 
these assumptions all do hold simultaneously, thus showing that the use of Black's model does 
not lead automatically to arbitrage opportunities. 

□ 



Example 6.3 Caps: Lognormal Bond Prices 

A cap can be decomposed into a portfolio of puts on zero coupon bonds. To be precise, the 
n th caplet (from the previous example) has 



Payoff = A5 n [L(T n _ 1 ) - R} + at time T, 



Since L(T n _i) is known at time T n _i this is equivalent to a time-T n _i payoff of 

M f+TuT~T = A t 1 - & + 5 nR)P{T n -uT n )] + 
1 + 6 n L{T n -i) 

+ 



= A(l + 5 n R) 



l + 5 n R 



- P(T n -i,T n ) 



This last line is easily seen to be the time-T n _i payoff of a portfolio of A{1 + <5„i?)-many put 
options with strike i + \ nR and expiry T„_i on underlying security P(t,T n ). If at time T n _i 
the caplet has the same payoff as a portfolio of puts on P(t,T n ), then, by the Law of One 
Price, the value of the caplet must have the same value as the portfolio of puts at any earlier 
time as well. 

Thus the t = 0-value of the n th caplet is 

Cn(0) = ^4(1 + 5 n R) x value of put option on P(t, T n ) with strike 



l + S n R 
and expiry T n _i 



This can be evaluated using the method of the first example of this subsection. 

□ 



Example 6.4 Swaptions: Lognormal Swap Rates 

Suppose we initiate, at time t, a pay-fixed interest rate swap starting at time T > t, with 
tenor structure T = To < T\ < ■ ■ ■ < Tjy on a notional amount A. This is known as a forward 
swap or deferred swap. Let 5 n = T n — T n _i, and recall that at T n pay-fixed receives 

A<5 n (L(T n _i) - St, T ) n = l,...,N 

where St t T is the T-forward swap rate at time t, and L(T n _i) is the spot LIBOR rate at time 
T n _i for the period [T n _i,T n ]. Further recall that St t T is the rate which sets the initial (i.e. 
time t) value of the forward swap equal to zero. 

The interest payments on a pay-fixed swap are equivalent to the payments of a portfolio 
consisting of short a coupon bond with coupon rate S^t, and long a floating rate note. The 
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bond and the FRN both come into existence at time T. The current value of such a forward 
starting bond bond is 

N 

A[J2SnSt, T P(t,T n ) + P(t,T N )} 

71=1 

The floating rate note will trade at par at time T, i.e. we need to set aside AP(t, T) at time 
t to purchase the FRN at time T. Hence the forward swap rate satisfies 

N 

-A[J2 S n P(t, T n )S t , T + P(t, T N )} + AP(t, T) = 

n=l 

(where the coupon bond and FRN have the same payment dates as the swap, and the same 
notional) and thus 

= P(t,T) - P(t,T N ) 

^ En=l6nP(t,T n ) 

If t = T, then St t t is just the ordinary spot swap rate at time t. 

A swaption C is the right to enter into a pay-fixed swap at some future date T at a strike 
rate R. If the tenor structure is T = Tq < T± < T2 < ■ ■ ■ < T/v, then the swaption gives the 
holder the right (but not the obligation) to receive at each of the dates T\ , . . . , T/v an amount 

A5 n (L(T n _x) - R) 

If a pay-fixed swap were to be entered at time T at the spot swap rate, then payments would 
be 

A6 n (L(T n -i) - St,t) 

and thus the swaption would be exercised only if R < St,t- The swaption thus gives rise to 
a series of payments 

AS n (S T , T - R) + 

at times T n . Each payment is equivalent to the payoff of A5 n -many calls with strike R 
and maturity T on underlying St,t- Using the generalized version of Black's model, i.e. 
assuming that St,t is lognormal under the riskneutral measure and making the appropriate 
approximations, the t = 0-value of each such payment is 



AS n P{0,T n ) [5 , T iV(di) - RN{al 2 )} 

<tVT 

St,t- Hence the value of the swaption is 

TV 



In S ° ,T I 1 a 2 T 

where d\ = R 7= , d 2 = d\ — a^T, and a is the volatility of the future spot swap rate 



C = Y, A5 n P(0, Tn^rNid!) - RN(d 2 )] where 



n=l 



S 0,T 1 l„2r 



In + ko l T 

1 



d = '' 

o-Vt 

d 2 = d 1 - aVf 



and 



So,T - 



P(t,T) - P(t,T N ) 
E!=iV(0,T n ) 
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We saw that we can make the Black formula for caps exact, provided we work with the 
appropriate numeraires, under the appropriate equivalent martingale measures. Can we make 
Black's formula for swaptions exact? Yes, indeed. Note that the numerator in the expression 
for So T is equivalent to a portfolio of zero coupon bonds, i.e. 



N 



^<5„P(0,T n ) 



n=l 

corresponds to a stream of cashflows of size 8 n at time T n . If, as is often the case, all the 8 n 
are of the same size, then this portfolio is just an annuity. Now we may think of the portfolio 
as a traded asset, call it X, and use it as numeraire. 

The first of the Black approximations is exact under the measure <Qx '■ The time-T value 
of all the payoffs of the swaption is 

N 

C T = Y, AS n P(T, T n )[S T , T -R} + = AX t [S t ,t - R} + 

71=1 



Hence 



so that 



Co 



Ct 
Xt 



Eq x [A(St,t-R)- 



N 



C = J2 M nP(Q,T n )E Qx [(S' 



T,T 



RY 



i.e. we discount outside the expectation. 

As for the second approximation, we need to show that the forward swap rate «So,t (which 
can now be seen to equal p ( > T )~ p (°> Tjv ) ) jg j us ^ the expected value of the future spot swap 



rate S T t under the EMM 



x 
\x, i.e. 



that Eq x [5t,t] = So t T- To see this, consider a pay- 



fixed forward swap F initiated at t = to start at time T, with interest payment dates 
Ti, . . . ,Tjv- The t = 0-value of the contract is Fq = 0, whereas at time T the value is 
Ft = Yln=i A5 n P{T,T n )[S T ,T - S ,t] = AX t [S t ,t - So,r]- The desired result now follows 
immediately from the fact that 



= |° = E Qx 
Xq 



F t 
Xt 



Hence Black's formula is exact, provided we assume that swap rates are lognormally dis- 
tributed under the EMM associated with the annuity process X t = J2n=i ^nP{t,T n ). 



□ 



It's pretty amazing that the Black formula for various derivatives (published in 1976) can 
in many cases be made exact using the change of numeraire technique (discovered in the early 
1990's). In particular, both the Black formula for caps and that for swaptions are exact if we 
assume that LIBOR rates are lognormal under the appropriate forward riskneutral measures, 
and that swap rates are lognormal under the "annuity" measure. 
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6.2 Review of Changes of Measure and Numeraire; LIBOR Rates 

Fix a horizon T* > and suppose that (fi, T, P, (<Ft)t, {S\)ij) is a market model, where the 
filtration (^t)t is generated by a standard (multi-dimensional) P-Brownian motion (Wf)t, 
augmented to satisfy the usual conditions. Let Q be the riskneutral measure, i.e. a measure 
which has the property that all asset price processes SI are martingales when denominated 
in units of the money market account A t . We briefly recall some facts about how Girsanov's 
Theorem is used to change the measure (e.g. to construct Q from P): 

• Assume that the asset dynamics are given by 

— * = p% S t ) dt + a% S t ) dW t —* = r t A t dt 
with suitable initial conditions. 

Recall that the market price of risk Af is a vector satisfying 

(This looks like it depends on the asset S l , but we know from previously developed 
theory that, for a model to be arbitrage-free, all assets must have the same market 
price of risk. Hence we've suppressed an index i.) 

• Let u(t, uj) be a predictable process, to be used as a kernel for a Girsanov transformation. 

• Define a new measure P by 




= e tf Ut d w t -yT \\u t \\ 2 dt 



• Girsanov's Theorem states that 

W t = W t - f u s ds 
Jo 

is a P-Brownian motion. 

• Thus the new asset dynamics are, under P, given by 

f iqi a A 

— * = (/4 + a\u t ) dt + a\ dW t -± = r t A t dt 

It follows that the market price of risk under P must satisfy the relation 

>h4 = ti + o\u t -r t = (Af + u t )a\ 

and thus 

Af = Ar + n, 

• Hence a Girsanov transformation adds the Girsanov kernel to the market price of risk. 
It adds volatility x kernel to the drift. 
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• To obtain a riskneutral measure Q, the new market price of risk Aj must be zero, and 
thus we must have ut = — Af . This is in agreement with what we found earlier. In that 
case, the drift becomes \x % — <r|Af = r t , which we already know very well. 

• To change from the riskneutral measure Q to an equivalent martingale measure Qx 
for numeraire X, we proceed as follows: Start in the riskneutral world, where = 

r dt + a s dW^, and ^ = r dt + a x dW^. Under Q x , the ratios S t = jfc are 
martingales. Now under Q, the ratios have dynamics 

^ = -vx(<?s ~ ox) dt + (a s - a x ) dW^ 
St 

= -a x o dt + a dW^ 

(where a = as— ax)- To make the drift equal to zero (i.e. to make St into a martingale), 
we need to to add a x ax = volatility xax, i-e. we need to use a Girsanov transformation 
with kernel ax- Thus 

• Hence we need to add ax to the riskneutral market price of risk to obtain the market 
price of risk under Qx- Since the riskneutral market price of risk is zero, the market 
price of risk under Qx is just the volatility of the numeraire X. 

• Numeraire-denominated asset price dynamics under the associated equivalent martin- 
gale measure are therefore just 

< f = (a s - a x ) dW^ 
ot 

• If the numeraire is the T-bond P(t,T), the associated EMM is called the T-forward 
riskneutral measure, and denoted by Q T . If bond price dynamics are 

^§=Vs(t)dt + a s (t) dWf 

under the "real-world" measure P, then the numeraire denominated dynamics are given 
by 

dP(t,S) T 

— = (f s - 0t) dW f 

P(t,S) 1 

where P(t, S) = p^'j^ and Wj is a Q T -Brownian motion. 

• Given future times T < S, the market price of risk under Q s is just as, whereas the 
market price of risk under Q T is ax- To move from Q s -world to Q T - world, we must 
change the market price of risk from as to cjt, i.e. we need to add ax — as to the market 
price of risk under Q s . Hence the change from Q^-world to Q T -world is effected by a 
Girsanov transformation with kernel ax — as, i.e. 

w =£t \Jo (JT ~ (JS * 
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We can also verify this directly. Recall that the Radon-Nikodym process & = Eqs [^p-|P<] 
for a change of numeraire is given by a ratio of asset ratios: 

_ P(t,T)/P(t,S) 
? * P(0,T)/P(0,S) 

and thus 

ddt = t;t[<TT(t)-*s(t)} dwf 

The solution of this SDE, together with the initial condition £o = 1> is just £ t = 
S t (j a T (u)-a s (u) dWf). 

Finally, note that the asset ratio process P(t, T) = satisfies the same SDE as does 

= - a s (t)] dwf 



■ - • - cru(l)\ f/ll'," ;,S 

P(t,T) 

although their initial conditions may differ. Hence £t and P(t, T) differ by a constant 
factor, i.e. 



St = cP(t,T) = 



cP(t,T) 
P(t,S) 



Let T* > be a a horizon for our bond market model. The time-t forward LIBOR rate 
for the future interval [T, T + 6} (where T < T* - 5) is defined by 

We saw earlier that L(t, T) is the interest rate for the period [T, T + 5] that can be locked in 
at time t (by a judicious investment in a portfolio of T- and T + <5-bonds with zero initial 
cost). 

Alternatively, the forward LIBOR rate L(t, T) can be regarded as the swap rate for a 
single-period swap settled in arrears. For suppose that we have a single-period interest rate 
swap, contracted at time t, for the period [T,T + 5], to be settled at time T + 5. Thus, at 
time T + 5, the pay-fixed side pays SR, and the receive-fixed party pays P _1 (T, T + 5) — 1, 
where P is the fair swap rate, and P _1 (T, T + <5) = 1 + 5S, S = L(T, T) the spot rate at time 
T for period T,T + 5]. Equivalently, by adding 1 to both payments, pay-fixed pays Y^ x and 
receive-fixed pays Y? 1 , where 

Y fx = l + SR Y fl = P _1 (T,T + 5) 

We can regard Y' x and Y* 1 as contingent claims which are paid out at time T + 5. It is clear 
that the time t-value of Y^ x is just 

Y fx = P(t,T + 6)[1 + 5R] 

The time-t value of Y? 1 is obtained as follows: If, at time T, we invest $1.00 in T + 5-bonds, 
the payoff at time T + 5 will be P _1 (T, T + 5). To obtain the required $1.00, we must invest 
in one T-bond at time t <T. Hence 

Y/ l = P(t,T) 
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The swap rate at time t is the rate R for which Y-f' x = Y-f l , and thus R = r+'^ +< ^ = 

L(t,T). 
Define 

P(t, T, S) = ffi ^ = 1 + <SL(t, T) for i < T < 5 and 5 = S - T 
P(t, S) 

Then P(t, T, S) is a Q s -martingale. In particular, the LIBOR forward rate L(t, T) is a Q T+5 - 
martingale. Thus the LIBOR forward rate L(t, T) is simply the expected value of the LIBOR 
spot rate L(T, T) at time T, where the expectation is taken under the (Q) T+<5 -measure. 



7 Lognormal Forward LIBOR Market Models 

We start with a pre-specified sequence of times 

= T < Ti < T 2 < ■ ■ ■ < T N = T* 

These times, typically settlement- or reset dates, are collectively known as the tenor structure. 
We also define Sj = Tj — Tj-i for j = 1, . . . , N. Then the forward LIBOR rate satisfy 

1 + S j+1 L(t,Tj) = JL ( *' Tj "\ = P(t,Tj,T j+1 ) 

We assume that the bond market satisfies a strong form of the no-arbitrage condition, i.e. 
we assume that there exists a riskneutral measure Q simultaneously for all discount bonds 
P(t,T). We denote, for each P(t,T), its associated forward riskneutral measure by Q T . Wt 
and Wj will denote, respectively, Q- and Q T -Brownian motions. 

Let S(t, T) be the volatility of the T-bond P(t, T) at time t. From the previous subsection, 
we know the following: 

• Q T i is obtained from Q T ^' +1 via a Girsanov transformation with kernel S(t, Tj)—S(t, Tj + i), 
i.e. 

= S T] ( S(t, Tj ) - S(t,T j+1 ) dW T t - 



• Each asset ratio P(t,Tj,Tj + ±) = is a Q T J +1 -martingale. 

• Each forward LIBOR rate L(t,Tj) is a Q T J +1 -martingale. 



The Q t j'+! -dynamics of the asset ratio P(t,Tj,Tj + i) are 



f(t,lj,lj + i) 



There is a constant c such that the Radon-Nikodym process and the asset ratio process 
are related 



3 + 1 



= cP(t, Tj,Tj + i) = c(l + 6 j+1 L(t, Tj)) 
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Note that, assuming that the forward LIBOR rate processes L(t,T) are strictly positive, we 
have the following dynamics: 



dL(t,Tj) = L(t,Tj)\(t,Tj) dW t 



This follows from the Martingale Representation Theorem: L(t,Tj) is a Q 7 J+ 1 -martingale, 
and thus we must have dL(t,Tj) = h t dW t J+ . Since L{t,Tj) is strictly positive, we may 
define X(t, Tj) = frach t L(t,Tj) to obtain dL(t,Tj) = L(t, Tj)X(t, Tj) dwj 3+1 . 



Now P(t,Tj,T j+ i) = l+S j+1 L(t,Tj), so that dP(t, Tj, T j+1 ) = S j+1 dL(t,Tj) = S j+1 L(t,Tj)X(t,Tj) dW t 



We also found that d p^x J T J ^) = ~~ S(t,Tj + \)) dwj^ 1 , and equating these expres- 

sions yields 

This expression will play an important role in the inductive construction of lognormal models 
of forward LIBOR rates. 

Since the move from (Q^^ 1 -world to Q^'-world is accomplished by a Girsanov transfor- 
mation with kernel S(t, Tj) - S(t, T j+1 ) = , the dynamics of L(t, Tj) under 
are given by 

dL(t, Tj) = L(t, Tj) \ ^ K )^ ] l\ dt + \(t, Tj) dW? 
[1 + d j+ iL{t,lj) 

because volatility x kernel must be added to the Q T J +1 -drift of L(t,Tj), while leaving the 
volatility unchanged (and the drift is zero, while the volatility is X(t,Tj)). 

7.1 The Brace— Gatarek— Musiela Approach to Forward LIBOR 

In most markets, caps and floors form the largest component of an av- 
erage swap derivatives book. ... Market practice is to price the option 
assuming that the underlying forward rate process is lognormally dis- 
tributed with zero drift. Consequently, the option price is given by the 
Black futures formula, discounted from the settlement data. 
In an arbitrage-free setting, forward rates over consecutive intervals 
are all related to one another, and cannot all be lognormal under one 
arbitrage-free measure. That is probably what led the academic com- 
munity to a degree of skepticism toward the market practice of pricing 
caps. . . 

The aim of this paper is to show that market practice can be made con- 
sistent with an arbitrage-free term structure model. . . This is possible 
because each rate is lognormal under the forward (to the settlement 
date) arbitrage-free measure rather than under one (spot) arbitrage- 
free measure. Lognormality under the appropriate forward and not spot 
arbitrage-free measure is needed to justify the Black futures formula with 
discount for caplet pricing. 

- Brace, Gatarek, Musiela [1997] 

The BGM-model starts from a family P{t, T) of discount bond prices up to some horizon 
maturity T* . We assume that each forward rate is over a period of length 5 (the same for all 
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rates). The bond price processes also give us the bond ratio processes (i.e. forward prices) 
P(t,T,S) = ^gj. The forward LIBOR rates L(t,T) are thus defined by 

1 + SL(t, T) = P(t, T,T + 5) for T < T* - 5 

BGM put their model inside the HJM framework, i.e. they assume that a term structure 
of instantaneous forward rates for all maturities (less than the horizon date T*) is available. 
In contrast, the Musiela-Rutkowski and Jamshidian approaches require forward rates only 
for a discrete set of tenor dates, as we shall see. Now recall that if, in an HJM model, the 
riskneutral dynamics of the instantaneous forward rate f(t, T) is given by 

df(t, T) = a(t, T) dt + a(t, T) dW t where a(t, T) = a(t, T) a(t, u) du 

(using the HJM drift condition), then the riskneutral bond price dynamics are given by 

d SrS =r t dt + S(t, T) dW t where S(t, T) = — f a(t, u) du 

Further recall that earlier we obtained 

6L(t,T) 



l + SL(t,T) 



X(t,T) = S(t,T) - S(t,T + 5) 



(which also follows if we apply Ito's formula to the identity 1 + 8L(t, T) = e^r + du and 
compare the dW^-terms). The main problem is this: 

How can we specify bond volatilities S(t, T) (or equiva- 
lently, the instantaneous forward rate volatilities a(t,T) = 
— 9S gj T ' > ) so that the resulting discrete simple forward LI- 
BOR rates will have the desired deterministic volatity struc- 
ture? 

We have already seen that L(t, T) is a non-negative Q^+^-martingale. For 1 + SL(t, T) = 
P(t, T,T + S), and so dL(t, T) = 5~ l dP(t, T,T + 5). But P(t, T, T + S) is a Q T+<5 -martingale 
(by definition of Q T+S ), with dynamics ^Itt+S) = [S(t,T)-S(t,T+5)} dW? +s . It therefore 
follows that 

dL(t, T) = 5- l P{t, T,T + S)[S(t, T) - S{t, T + 5)} dWj +h 

= L( t ,T)(' 1 + ^ r "^fj-^ T + { ") mr> 

i.e. 

dL(t,T) = L(t,T)X(t,T) dWj +5 

where A((>r) _ [1± mxm^hl±3 

We are therefore able to derive the forward LIBOR dynamics directly from the bond price 
volatilities (or, equivalently, the instantaneous forward rate volatilities). Since the forward 
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riskneutral measure Q T+S is obtained from the (spot) riskneutral measure Q by a Girsanov 
transformation with kernel S(t, T + 5), we have 

dW? +s = dW t - S(t, T + 5)dt 

for a Q-Brownian motion Wt- Thus the riskneutral drift is directly determined by the volatil- 
ity structure (as it is in the HJM model), giving riskneutral forward LIBOR rate dynamics 

= -X(t, T) ■ S(t, T + 5)dt + X(t, T) dW t 

Now suppose that we want to create an HJM model in which forward LIBOR rates L(t, T) 
have a deterministic volatility structure X(t,T). Above, we found that 

S(t, T) - S(t, T + 5) = £ +S a(t, u) du = _i^L_ A(t, T) 

(where S and a are the bond and instantaneous forward rate volatilities respectively). In 
order to find the bond volatilities, it is necessary to impose some additional conditions. Set 

a(t, u) = when < u — t < 5 

(This is the fundamental assumption made in BGM(1997)). 
Now find the bond volatilities by a recursive procedure: 

• Choose n such that n5 < T — t < [n + 1)5. Equivalently n = sup{k £ N : kS < 
T — t} = [(5 _1 (T — t)] (where [x] is the integer part of x). 

• Then S(t, T — n8) = — f^~ n& cr(t, u) du = 0, because < u — t < 5 when t < u < 
T-nS. 



Thus 



S(t, T) = [S(t, T) - S(t, T-5)} + [S(t, T-S)- S(t, T - 26)} + ... 

... + [S{t, T-{n- 1)5) - S(t, T - n5)} 



implies 



o/j. rp\ 5L(t,T-5) 5L(t,T-25) 

S{t > T) = -l + 5L(t,T-5) Kt > l + 5L(t,T-25) Kt > T ~ 26) ~ " ' ' 

SL(t,T-nS) w m 

-X{t,T-n5) 



l + SL(t,T-n5)' 
• i.e. 

k=l 

Equivalently, 
(i) Define S(t, T) = for < T - t < 5. 



Interest Rate Modelling 



65 



(ii) Then define S(t, T) = S(t, T - 5) - ^I'^t-S) ^ T - 5) ioi 5 < T - t < 25. 

(Note that if <5 < T - t < 28, then < (T — 5) - t < 5, so S(t, T — 8) has already been 
defined.) 

(iii) Then define S(t, T) = S(t, T - 5) - ij§ftSS) A(t, T - 5) for 25 < T - t < 35. 

(Note that if 25 < T - t < 35, then 5 < (T - 5) - t < 25, so S(t, T - 5) has already been 
defined.) 

(iv) . . . etc. 

In this way, if we specify bond volatilities by this forward induction, then we will have an HJM 
model in which the forward LIBOR rates L{t, T) have the required deterministic volatilities 
\{t,T). Since each L(t,T) is a strictly postive Q T+<5 -martingale, it follows that each L(t,T) 
is lognormal under Q T+<5 , and thus that the Black formula for caps is valid in this model. 

7.2 The Musiela-Rutkowski Approach to Forward LIBOR 

Unlike the BGM-approach, which lies within the HJM framework and specifies a model 
of forward LIBOR rates L(t,T) for all maturities T (below the horizon T*), the Musiela 
Rutkowski (MR) approach only specifies LIBOR rates for a discrete set of maturities. We 
start with a discrete tenor structure 

< T < Ti < • • • < T N = T* 5 n = T n - T n _! 

and define T_i = (for ease of handling certain formulas). We further assume that we are 
given 

• A family of bounded adapted processes X(t, T n ) for n = 0, . . . , N — 1 which represent 
the volatilities of the forward LIBOR rates L(t,T n ). 

• An initial term structure P(0, T n ) of discount bond prices (used to specify the initial 
conditions of the SDE's which we will write down for the LIBOR rates). We further 
assume that P(0, T ) > P(0, T x ) • • • > P(0, T N ). 

In contrast to the BGM approach, we do not need a bond price dynamics at all, i.e. we will 
attempt to model LIBOR rates directly. 

Before we construct the MR model of LIBOR rates, a lemma which will prove useful 

Lemma 7.1 If X,Y are adapted processes 

dX t = a t dW t dY t = (3 t dW t 

and if Z t = i^y t , then 

d(Z t X t ) = Z t (a t - foZ t X t ) ■ (dW t - foZ t dt) 
i.e. d(Z t X t ) = r) t ■ (dW t - foZ t dt) 

for some process r/t. 

Proof: A straightforward application of Ito's formula. 
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□ 

Whereas the BGM approach shows how to define bond volatilities by forward induction, 
the MR approach directly constructs a set of measures under which forward LIBOR rates 
have the required volatility structure by backward induction. It is therefore convenient to 
introduce the following backward notation. Put 

T fc * = T N ^ k so that T* = T* > T* > ■■■ > = T 

We start by working under a T^-forward riskneutral measure Q Tjv = Q T( *, together with 
a Q TiV -Brownian motion W Tn = W t . it is not necessary to construct this measure: we 
can assume that Q Tn is the measure P which governs our model, and that W Tn is the Wt 
which drives the economy Ultimately, we will be able to specify all the dynamics under this 
measure, the terminal measure. Let L(t,T 1 *) = L(t,T^-i) be a process which satisfies the 
SDE plus initial value 

dL(t,T?) = L(t,T^)X(t,T[) dWf N 

T(n T », P(o,?T)-p(o,r *) 

This defines the forward LIBOR rate L(t, T*) = L(t,T N -i) in the MR model. 

We now use this to define the forward LIBOR rate L(t,T£) = L(t,Tjv-2)- To do so, we need 

to construct the forward riskneutral measure for maturity T| . Under Q T 2 5 a ll the bond ratios 

Pit T* ) 

p| t ' r ") are martingales. Now define the ratio 

Pit Th) Pit T*) 

U N - n+1 (t, T k ) = - or, equivalents U n (t, T fc *) - ' ' 



P(t,T n ) ' ^ J ^ P^T* n _ x ) 

and note that each U n (t,T£) is required to be a martingale under the measure Q 7 ™- 1 (which 
we must still construct). Further note that 



l + S N L(t,T*) 

so that by the lemma, 

for some process r/k,t (whose exact nature is not important right now). In order for each 
U2(t, T£) to be a martingale, it suffices to find a measure under which 

is a Brownian motion. This is possible if we perform a Girsanov transformation from Q Tjv = 
Q T o with kernel 7 (s, T*) = ^^g^ A(s, T*), i.e. if we define 



dQ T o 
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We now let L(t,T|) be a process which solves the SDE and initial condition 

dL(t,T%) = L(t,TZ)\(t,T%) dW T t ' 

We continue in this way: Suppose that we have already constructed the LIBOR rate processes 
L(t,Tf), . . . ,L(t,T*), for n < N — 1. Suppose further that this has been done so that each 
forward measure and Brownian motion has been specified, in particular that we have already 
constructed Q T "-i and and that dL(t,T*) = L(t,T*)X(t,T*) dW^- 1 under Q T »-i. 

We must now construct a measure Q T ™ and an associated Brownian motion wj n . We require 
that each U n +i(t,T£) is a Q T ™ -martingale. Now 

Un+l[t ' Tk) -l + 5 N - n+1 L(t,T*) 
Using the lemma, we see that 

du n+1 ( t ,m = m ■ (dw T ^ - ^p^I^L k^t:) dt ) 

\ 1 + d N - n+ iL(t,l*) J 

for some process 77^ (whose exact nature is not important right now). In order for each 
U n +i(t,T%) to be a martingale, it suffices to find a measure under which 

1 1 Jo 1 + S N - n+1 L(s,T*) v n) 

is a Brownian motion. This is possible if we perform a Girsanov transformation from Q T ™-! 
with kernel j(s,T*) = rfe^ggj|) A(s, T*), i.e. if we define 

= S T * ( [ i(s, T*) dW?"- 1 

We now let L(t,T* +l ) be a process which solves the SDE and initial condition 

dL(t,T* +1 ) = L(t,T* +1 )X(t,T* +1 ) dwP 

naT * ,_P(Q,r* +1 )-p(o,r*) 

( ' n+1 '~ <W(0,T*) 

We have now constructed a sequence of processes L(t, T n ) which are models of the forward 
LIBOR rates, with the desired volatilities. Since we also know the Girsanov kernels of each 
transformation, we can specify all LIBOR rate dynamics under the terminal measure. Induc- 
tively, 

dL(t, T*) = L(t,T*)X(t,T*) dwT^ 1 

= —L(t, T*)X(t, T*)~/(t, T*_ ± ) dt + L(t,T*)X(t,T*) dW^ 2 

= -L(t,r*)A(t,r*)[ 7 (i,T„*_ 1 ) + 7 (t,r*_ 2 )] dt + L{t,T*)x{t,T* n ) dw^' 3 

n—l 

= -L(t,T*)X(t,T:)J2l(t,T*_ k ) dt + L{t,T*)X{t,T* n ) dW? 
k=i 
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where 



rri* \ _ 5 N _ k+1 L(t,T£) 

7(t ' Tfc) -i + ^- fc+ iL(t ; r*) A(t ' Tfc) 



and hence, when we translate from backwards time to ordinary time, 

The Musiela-Rutkowski forward LIBOR rate dynamics under the ter- 
minal measure Q Tn are given by 

dL(t,T n ) = -L(t,T n )X(t,T n ) £ 5 ^f T ^ T T ^ dt+L(t,T n )X(t,T n ) dW?» 

This must be solved recursively: First find the solution for L(i,T/v-i). 
Once this has been found, find the solution for L(t, Tiv-2)- Note that the 
SDE for L(t, Tjv-2) also contains L(t, Tjv-i), but we've already found 
that. Then solve the SDE for L(i,T/v-3) (which contains L(t, T n _i) and 
L(t, Tjv-2); these have been determined). And so on. . . 

It is therefore possible to find a model in which LIBOR rates have the required volatilities 
\(t,T n ). If these volatilities are deterministic, then each L(t,T n ) will be lognormal under 
Q Tn+1 . In that case, the Black formula for caps will be exact. 

7.3 Jamshidian's Approach to Forward LIBOR 

Like the Musiela-Rutkowski approach, Jamshidian(1997) does not require bond price dynam- 
ics, and models LIBOR rates for a discrete set of tenor dates = T_i < Tq < T\ < ■ ■ ■ < 
T/v = T* via a backward induction. But instead of working under the terminal measure, 
Jamshidian defines a spot LIBOR measure. This measure is obtained if we take as numeraire 
a certain portfolio of zero coupon bonds with unit initial value. 

We begin by observing that the prices of discount bonds are not completely determined 
by the forward LIBOR rates. This is true at tenor dates, but if t lies between tenor dates, 
e.g. T n < t < T n+1 , then P(t, T n+k ) = P(t, T n+1 ) ■ 1+5n+2 l {t>Tn+l) i+WmW^) - Thus 
knowledge of the LIBOR rates is not enough — we also have to know the discount factor 
to the next tenor date (i.e. P(t,T n+ i)). By working under the spot LIBOR measure, this 
problem can be circumvented. 

Consider the following portfolio of discount bonds X. Its initial value is $1.00. At all 
subsequent times, all wealth is invested in the next-to-mature bond. Thus at t = 0, $1.00 is 
invested in P(t,To). At To, the payoff of these bonds is reinvested in P(t,T\) and at Ti, the 
payoff is reinvested in P(t, T2), etc. Thus at time T n , the value of the portfolio is 

P(T n ,T n+ i) 

X T n = 



P(0, T ) • P(T , Ti) P(T n , P(T n+1 ) 

= value of T„ + i-bonds x no. of T n+ i-bonds 



An instant later, when T n <t < T n+ i, the value X t of the portfolio is simply 

P(t,T n+1 ) 



X t = 



P(0, Tb) • P(T , Ti) P(T n , P(T n+ i) 
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because the value of the T„ + i-bond has changed, but the number of T n+ i-bonds in the 
portfolio has not. Hence 

n(t) 

x t = p(t,T n{t) ) ■ n p-^T^n) 

where nit) = inf {n : T n > t} 

A spot LIBOR measure Qx is obtained by taking X t as numeraire, so that each asset ratio 
process p ^ n ^ is a Qx-martingale. The asset ratios can be written as 

P(t,T n+1 ) _ P(t,T n{t) m n k=n{t)+1 (l + 5 k L{t,T k ^))-i 



n(t) n 

[JfiiWu.rn))" 1 n ^ + W,T k _ 1 ))- 1 

k=0 fc=n(t)+l 



n 



+ S k L(t AT^Tk^))- 1 



k=0 

Hence the prices of the asset ratios are completely determined by the LIBOR processes. 
We now aim to describe the LIBOR rate dynamics under the spot LIBOR measure Qx, 
and that this requires knowledge only of the LIBOR rate volatilities (and not, say, bond or 
instantaneous forward rate volatilities as well). For the moment, assume that bond price 
dynamics are given by some ltd processes 

d pQTj = m ^ Tn) dt + S ^ Tn) dWt 
under the "real-world" probability measure P. By definition of Xt (i.e. by (*)), we see that 

^ = m(t, T n(t) ) dt + S(t, T n(t) ) dW t 

Pit T ) 

Moreover, if we apply Ito's formula to 1 + 6 n+ iL(t, T n ) = p ^ , we see that 

m(t, T n ) - m(t, T n+l ) - (S(t, T n ) - S{t, T n+l ))S{t, T n+1 )) dt 



dL(t,T n )= PiLT " 



$n+lP(t, T n+ i 

- (s(t,T n ) - S(t,T n+1 j) dW t 
= fi(t,T n ) dt + ((t,T n ) dW t 

where 

fi(t,T n ) = - P ^\ T ^ (m(t,T n )-m(t,T n+1 ) ) - {(t, T n )S(t, T n+1 ) 

C(t,T n ) = P ^ (s(t,T n ) - S(t,T n+1 )) 
It follows that 

6 k+1 ((t,T k ) 



k=n(t) 
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for j > n{t). 

If 7 (i) is the Girsanov kernel for transforming P to Qx, i-e. if 

Pit T ) 

then r ^ t n > has zero drift under Qx- But 



£t* (J It dW t ) , 



d 



P(t,T n ) P(t,T n ) 



X, 



x. 



(m(t,T n ) - m(t,T n(t) ) - S(t,T n(t) ) ■ (S(t,T n ) - S(t,T n{t 
+ (s(t,T n ) - S(t,T n{t) )) dW t 



dt 



Now in the Girsanov transformation, (S(t,T n ) — S(t,T n ^)) ■ 7i is added to the P-drift to 
obtain the Qx^drift, which is zero, and so 

m(t,T n ) - m(t,T n{t) ) - S(t,T n{t) ) ■ (S(t,T n ) - S(t,T n(t) )) + {S(t,T n ) - S(t,T n{t) )) ■ 7t = 

which yields 

m(t, T n ) - m(t, T n+1 ) = (s(t, T n(t) ) - 7t ) ■ (s(t, T n ) - S(t, T n(t) )) 
for n = 0, . . . , iV. It follows that 

m(t,T n ) - m(t,T n+1 ) = (m(t,T n ) - m(t,T n(t) ) - (m(t, T n+1 ) - m(i, T n(t) 
= (s(t,T n(t) ) - 7t ) • (s(t,T n ) - S-(t,T n+1 )) 

Now multiply both sides of this equation by g +1 p( t T +1 ) ^° obtain 

5 n+ ^P^T n+1 ) H*' Tn) " T " +l) ) = C( *' Tn) ( S{t > Tn{t)) " 7 *) 
Looking back to the definitions of fi and £ in the dynamics of L(t, T n ), we see that 

H(t, T n ) = C(t, T n ) (S(t, T n(t) ) - 7t - T n+1 )) 

and hence 



dL(t, T n ) = C(t, T n ) 



(S(t, T n{t) ) - S(t, T n+l ) - 7t ) ^ + dW t 



These are, of course, the P-dynamics. To get the Qx _ dynamics, we must add volatility x 
kernel = Q • 7 t to the drift to obtain 



dL{t,T n ) = C(t,T n ) 



S(t,T n(t) ) - 5(t,T n+ i)) dt + dWj 



where = Wt — JqJu du is a Qx^Brownian motion. Finally, using (**), we obtain 



dL(t,T n )= £ 

fc=n(t) 



<y fc+ iC(t,T fc ).c(t,T n ) 



These are the forward LIBOR rate dynamics under the spot LIBOR measure. 
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8 Appendix: Girsanov's Theorem 
8.1 Motivation 

When pricing contingent claims, we use risk-neutral valuation: The t = price of a claim X 
is the risk-neutral expectation of its discounted payoff. 

X = E Q [X] 

The measure Q is not the same as the "real-world" measure P — we have to change the 
probability measure. 

Suppose that we start with real-world asset dynamics, e.g. a GBM 

^ = n dt + a dW t 
St 

on a filtered space (0, JF, P, F), where Wt is a (F, P)-BM. There are two questions that concern 
us: 

• What happens to the dynamics of St when we change measures? 

• How do we actually go about changing measures? 

Example 8.1 Earlier, we introduced Wiener space (C[0, oo), C, P), which is a probability 
space with sample space 0, = C[0,oo), the set of all continuous functions from [0, oo) to R. 
This space comes equipped with a stochastic process X, the coordinate process, defined by 
X t (u) = uj{t). (Remember that each u> £ C[0, oo) is a function.) The cr-algebra C is generated 
by this process. Wiener space also comes equipped with a measure P — Wiener measure - 
which has the property that the coordinate process X is a standard Brownian motion under 
P. 

Suppose now that we change the measure on C[0,oo), as follows: Let oj$ : [0,oo) — ► R : 
t i ^ be the constant function with value 0. Let Q = 8 Wo be the Dirac point mass, so that 
Q(A) = 1 if ujq G A, and Q(A) = otherwise. Now what are the dynamics of the coordinate 
process XI 

Clearly, Q({co : \/t(u>(t) = 0)}) = Q({o;o}) = 1, and hence, under Q, the process X is a.s. 
constant with value zero. 

So though X is a Brownian motion under P, it doesn't remotely resemble a Brownian 
motion under Q. 

Note however, that P, Q are not equivalent measures. It turns out that this sort of thing 
cannot happen when the measures are equivalent. 

□ 

Now recall that the stochastic integral 6t dSt has a very important financial inter- 
pretation: It is the gain made by trading the portfolio 6. Stochastic integrals can only be 
defined for integrators S which are semimartingales, i.e. it can only be defined if we can find 
a decomposition 

St = S + M t + A t 

where M is a local martingale, and A a finite variation process. It is rather clear that 
changing measures usually destroys the martingale property — think of a game of coin tossing, 
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first under a measure where the coin is fair, and then change the measure to one where 
P(Heads) > 0.5. So the decomposition above depends on the measure: The P-local martingale 
may not be a Q-local martingale. As a consequence, it is not obvious that a P-semimartingale 
is necessarily a Q-semimartingale, with the result that gains may be undefined in the risk- 
neutral world. 

It turns out that, provided P and Q are equivalent measures, all remains well. We now 
show why. 

8.2 Girsanov's Theorem — General Statement 

The following is adapted from Stochastic Calculus: A Practical Introduction by Rick Durrett, 
and Stochastic Integration and Differential Equations, by Philip Protter. 

Because we have defined stochastic integrals only for continuous semimartingales, we 
will assume that our filtration only admits continuous martingales (up to a modification). 
This is the case, e.g., for the filtration generated by Brownian motion (by the Martingale 
Representation Theorem). 



This is a good time to recall Bayes ' Theorem for calculating conditional expectations when 
we change the measure: If Q << P, £ = ^ and £t = Ep[£|jFt], then 



We say that two measures P, Q on a measurable space with filtration F = {Ft)t are locally 
equivalent iff P, Q agree on F t , for every t. Let F t ,Q t be the restrictions of P,Q to Tt, 
respectively, so that F t (A) = F(A) for A € Tt, that 



Let at be a continuous version of the Radon-Nikodym process Note that a is a 

(F, P)-martingale: For if A G T s , then Q S (A) = Q t (A). Now Q S (A) = f A a s dF s = f A a s dF t , 
and Qt(A) = f A at dFt, which implies that a s is an jF s -measurable variable with the property 
that 



and thus that Ep t [at |.F S ] = a s , by definition of conditional expectation. 

Lemma 8.2 A process Y t is a Q-(local) martingale iff a t Y t is a F-(local) martingale. 

Proof: Fix s < t. Then by Bayes' Theorem, 



This proves the result for martingales. Now use localization to obtain the result for local 
martingales. 



E F [Z\F t }=E Ft [Z\Ft] 



etc. 





Y is Q-martingale 
& a s Y s = a s E Qs [Y t \F s } = E ¥t [a t Y t \F s } 
44> aY is a P-martingale 
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Exercises 8.3 Fill in the details in the above proof. 

□ 

Theorem 8.4 (Girsanov's Theorem) If Y is a F '-local martingale, then 

Y = Y - a" 1 • [a,Y] 

is a Q-local martingale. 

Proof: Let At = aJ 1 d[a,Y] s , and assume, for the moment, that this integral exists. We 
want Y t — A t to be a Q-local martingale, so it suffices to show that a t (Y t — A t ) is a P-local 
martingale. By Ito's formula, using [a, A] = (because A is of continuous and of bounded 
variation), we have 

d(a t (Y t - A t )) = a t {dY t - dA t ) + (Y t - A t ) da t + d[a, Y] t 
= a t dY t + (Y t - A t ) da t 

Since the integrators Y, a are P-local martingales, so are the preceding stochastic integrals. 
Hence at(Yt ~ At) is a P-local martingale, as required. 

We omit the (technical) proof that At is well-defined for all t. It may be found in Durrett. 

H 

Corollary 8.5 If Y is a P-semimartingale, it is also a Q-semimartingale. 

□ 

The following result is very important for mathematical finance. 

Theorem 8.6 Let Y be a semimartingale (under and let H be predictable. The 

quadratic variations [Y] are the same under P, Q. The stochastic integrals (H • Y) are the 
same under P, Q 

Proof: We present only a heuristic outline of the proof: Recall that if P ~ Q, then con- 
vergence in P-probability is equivalent to convergence in Q-probability. Since the quadratic 
variation of Y is a limit in probability of sums of the form YljO ! t j+1 ~ Ytj) 2 , the quadratic 
variations remain the same under P,Q. 

As for the stochastic integrals: It is clear that if A is a finite variation process, then H • A 
is the same under both measures, as it is defined pathwise (i.e. u-hy-uS) as a Stieltjes integral. 
It therefore remains to consider H • M, where M is a P-local martingale. Any continuous 
local martingale can be made bounded, by stopping. In that case (H • M)t is defined as an 
I? limit of (H n • M)t, where the H n are simple processes. Now the simple integrals H n • M 
are defined pathwise, and therefore coincide under both measures. Equivalence of measures 
ensures that the L 2 -limits of the H n • M are the same under both measures. 
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8.3 Exponentials 

Recall the form of the Doleans exponential: 

S(M) t = e Mt ~^ 

We know that 

d£(M) t = £(M)t dM t 

and hence £(M)t is a local martingale whenever Mt is a local martingale. Doleans exponentials 
play a very important role in changes of measure: Recall that if Q is equivalent to P, then 
at = ^ is a strictly positive P-martingale. Suppose now that Dt is an arbitrary continuous 
P-local martingale. Define 

M t = In D + I DJ 1 dD s 
Jo 

Then M is a continuous local martingale also. Moreover, by Ito's formula, 

din A = A" 1 dD t ~ \ Dt 2 d[D] t = dM t - 1 - d[M] t 
from which it follows that 

D t = £(M)t 

In fact, we have almost shown: 

Theorem 8.7 If D is a strictly positive continuous local martingale, there exists a unique 
continuous local martingale M such that D = £(M). 

Proof: Existence has just been dealt with, so only uniqueness remains to be shown. If 
M,N satisfy £(M) = D = £(N), then taking logarithms yields M - N = \([M] - [TV]). 
The lefthand side is a continuous local martingale, and the righthand side an FV process. It 
follows that M — N is constant (as the only continuous local martingales of finite variation 
are the constants), and then that M = N. 

H 

Example 8.8 Consider a one-stock Black Scholes model, with P-dynamics 

dS t = S t [(n-r) dt + adW t } 

over a finite horizon [0, T]. We want to find a measure Q under which St is a martingale for 
t < T. The way to do this is to find a continuous local martingale at for which at§t is a 
P-martingale. If we can do that, we may define Q on Tt by ^jjl = cut- 

So let's look for such an a: Noting that a$ = 1, we can, by the martingale representation 
theorem, represent a as a stochastic integral: 

a t = 1 + / H s dW s 
Jo 

Now define K t = so that a t = 1 + J * a s K s dW s . Then dat = atK t dWt, and so 
a t = £(K*W) t . 
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It follows that 



d(a t S t ) = a t S t K t dW t + a t S t [((i -r)dt + a dW t ] + a t K t S t a dt 
= a t S t [(K t + a)dW t + ( f i-r + K t a) dt] 

So atSt will be a P-local martingale provided that K t = r —^ L - 
Thus we may define a risk-neutral measure Q on Tt by 

— = £[ •W) onf r 

dP V o~ Jt 

The only problem that remains is: Is it true that Q a probability measure, i.e. is Q(fi) = 
E ¥ [£(K»W) T ] = 1? 

Since £(K • W) is always a non-negative local martingale, it is always a supermartingale. 
Hence the above question can be rephrased as: Is it true that £{K • W) is a martingale? 

□ 

Exercises 8.9 Suppose that M is supermartingale. Show that it is a martingale iff EX t = 
EX for all t. 

□ 

As we have seen, a question that will be important in mathematical finance is the following: 

Given that M is a continuous local martingale, when is £ (M) a genuine martingale? 

As a matter of policy, we have often glossed over the technical differences between martingales 
and local martingales. Here, therefore, we will simply state two criteria that partially answer 
this question. See Protter for proofs. 

Theorem 8.10 (Kazamaki's criterion) Suppose that M is a continuous local martingale with 
the property that 

supE[e5 MT ] < oo where the sup is over all bounded stopping times 

T 

Then £(M) is a UI martingale. 

Theorem 8.11 (Novikov's criterion) Let M be a continuous local martingale, and assume 
that 

E[e^ M H < oo 

Then £(M) is a UI martingale. 
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8.4 Girsanov's Theorem for Brownian Motion 

Theorem 8.12 (Girsanov's Theorem for Brownian Motion) Suppose an n-dimensional pro- 
cess Y has P -dynamics 

dY t = ntdt + a t dW t (t < T) 

where W is a standard d-dimensional P-Brownian motion, fit(uj) G ~M. n ,at{uj) € M nxa! . Let 
Xt(oj) G M d be predictable. Define a measure Q on Tt by 

d ® r(\ urt 
Assume that Novikov's condition holds: 

e r e i/ ( rii A «ii 2ds i < oo 

Then: 

(i) Q is a probability measure on Tt- 
(ii) W t = W t - £ X s ds is a Q-Brownian motion. 
(Hi) The Q-dynamics ofY are given by 

dY t = {^t + <rth) dt + a t dW t 

Proof: (i) We have <Q)(fi) = E ¥ [£(X • W) T ] = £(X • W) = 1, by Novikov's criterion. Hence 
Q is a probability measure. 

(ii) Note that at = ^ t = £{X»W)t has dott = atXt dWt, so that, by Girsanov's Theorem, 

W t = W t - f a' 1 d[a, W] s = W t - f X s ds 
Jo Jo 

is a Q-local martingale. Now since f' Xt dt is a continuous FV-process, we see that 

[W i ,Wi] t = [W\W% = 5 ij t 

Hence W is a continuous Q-local martingale with the same covariance process as a d- 
dimensional Brownian motion. By Levy's characterization, W is a Q-BM. 

(iii) follows from the fact that dYt = fit dt + at dWt and dWt = dWt + Xt dt. 

H 

Remarks 8.13 This has important consequences: Suppose, under P, we start with "real 
world" GBM dynamics 

dS t = S t [fi dt + a dW t ] i.e. dS t = S t [(fi - r) dt + a dW t ] 

Suppose we now construct a new measure Q as above. This Girsanov transformation with 
kernel X adds aX to the drift of S, but does not change the volatility: 

dS t = S t [(fi-r + aX) dt + a dW t ] 
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For Q to be a risk-neutral measure, St must be driftless, i.e. 

a A = r — fi 

With only one asset, this translates to A = i.e. the Girsanov kernel is minus the market 

price of risk. 

The fact that a Girsanov transformation does not affect the volatility is also important: 
It implies that we can use real-world observations to estimate risk-neutral world volatility. 

□ 

9 Appendix: Correlated Brownian Motions 

When many assets are available in the economy, it is unrealistic to assume that these are 
all driven by only one source of noise. It would be equally unrealistic, however, to assume 
that all are driven by separate, indpendent, Brownian motions. Thus it becomes necessary 
to generate multiple correlated Brownian motions. 

Let's first consider a simpler case, where we are trying to generate not correlated Brownian 
motions, just correlated normal random variables, i.e. suppose that we want to generate mean 
zero normal random variables X±, . . . , X n with a specific covariance matrix E = (cXy). Here 
uij = Cov(Xi, Xj). 

You can check the following simple 

Fact: If (Xi , . . . , X n ) is a random vector with covariance matrix E and if A is an n x n-matrix, 
then the random vector 









c 


H 











has covariance matrix AY<A tr . 

Indeed, Cov(Y) = E[YY* r ] = ^E[XX* r ]A* r . 

□ 

Covariance matrices are necessarily symmetric positive semidefinite, and it is known that 
symmetric positive semidefinite matrices have a Cholesky decomposition, which means that it 
is possible to find a (real) lower triangular matrix C such that 

E = CC tr 

Note that if C is an arbitrary matrix, the CC tr is necessarily symmetric (obvious), and 
positive semidefinite: If x is a column vector, then x* r C is a row vector, with length given 
by ||x* r C|| 2 = (rc* r C)(x* r C)' r = ^ r CC tr x. Since the length of a vector is necessarily non- 
negative, CC tr is positive semidefinite. 

Thus any matrix A that can be written as A = CC tr is necessarily symmetric positive 
semidefinite. By the Cholesky decomposition, the reverse is also true. Indeed we can find a 
lower triangular C which does the trick. There is no deep mathematics behind this — we 



78 



Correlated Brownian Motions 



merely need to solve 

/ an an ■■■ a in \ I en 
ai\ a22 ■ ■ ■ a-2n &i\ C22 



\ 



(cu 



C21 
C22 



c„i\ 

C„2 



\a. 



a n 2 



\Cnl C„2 



J \o 



This system is easily solved: c\ x = an gives us c\\. c\\Cq.\ = 012 now gives us C21, etc. 

There are fast algorithms available for calculating Cholesky decompositions. 

Now suppose we are able to generate independent standard normal random variables 
Zi, . . . , Z n . These have the identity matrix as covariance matrix. Define a random vector 
X = CZ. Then the covariance matrix of Z is simply CIC tr = CC tr = E. Thus to get a vector 
X of mean zero normally distributed random variables with covariance matrix E, proceed as 
follows: 

• Generate a vector Z of independent standard normal random variables (of the same 
dimension as X). 

• Find the Cholesky decomposition E = CC tr of the symmetric positive semidefmite 
matrix S. 

• Put X = CZ 

Note that if E = (cr^ ) is a covariance matrix, then the correlation matrix is given by 

Pij 



33 



Clearly the correlation matrix is also symmetric. 

Now to obtain correlated Brownian motions W l lt, ■ ■ ■ , W^, we can proceed in a similar 
way. But first: What exactly do we mean if we say two Brownian motions W 1 ,W 2 are 
correlated? Clearly this has meaning if we speak about changes in the processes. If W 1 , W 2 
are highly correlated, then we expect a positive change in W 1 to be accompanied by a positive 
change in W 2 . 

Now suppose that we have independent standard Brownian motions B\^B^. Consider a 
matrix V = (7^) with the property that all the rows of V have unit length. Define 



( w t\ 








= r 











so that each W\ = Y^jlijBt ls a linear combination of Bps. It follows that each W\ is a 
continuous local martingale. Now 



E^j (e^) t^'^ 
E^ 
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because [B^,B k ] t = Sj^t and Yltjlij = 1- Hence, by Levy's Characterization, each W\ is a 
Brownian motion. Now 

t = ^2 likljtfklt 



which we may also write as 
Hence 



dw l t dwi = (rr' r )y dt 



e[w?w/] = e[ [w\ w%] = (rr*%t 

Thus W l , W 3 are correlated Brownian motions, and the correlation between W\ and W/ is 
simply (rr tr )jj, independent of t (because the variance of each W\ is just t). 

Note that if X, p are, respectively, the covariance and correlation matrix of (W/, . . . , W t n ), 
then £ = pt. Hence p is also symmetric positive semidefinite, and thus has a Cholesky 
decomposition p = TT tr . 

Further note that not every symmetric positive semidefinite matrix can be the correlation 
matrix of some multidimensional Brownian motion: Since the correlation of a random variable 
with itself is 1, it is necessary that a correlation matrix has l's down the diagonal. This, in 
turn, implies that the Cholesky decomposition matrix T will have row vectors of unit length. 

Hence, to create correlated Brownian motions with correlation matrix p, proceed as fol- 
lows: 

• Find the Cholesky decomposition p = TV . T will have rows of unit length. 

• Define W = TB, where B is a multidimensional standard Brownian motion (with 
independent component processes). W will be a multidimensional Brownian motion 
with correlation matrix p. 

One final remark about differential notation: Since rr* r = p, and since dW\ dW 3 t = 
(TT tr )ij dt, we have 

d[W\ W j ] t = dW l t dWf = Pij dt 

Example 9.1 To create two correlated Brownian motions W^,W^ with correlation p (a 
number, not a matrix), proceed as follows: The correlation matrix is 

1 P 
P 1 

Its Cholesky decomposition is found by solving 

1 p\ _ fa 0\ fa b 
p \)~ \ b c) \0 c 



(= rr tr ) 



for a, b, c. (Recall that T is lower triangular.) Thus a = 1,6 = p,c = \J\ — p 2 . Gratifyingly, 
the rows of V are seen to possess unit length. 

Finally, if B},B 2 are standard independent Brownian motions, then 



W} = B\ 



Wf = P B\ + Vl - P 2 B 2 t 
are Brownian motions with correlation p. 
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Example 9.2 Suppose we have asset dynamics 



dS}\ _ f0.3S}\ t fO.lS} 0ASj\ (dWl 



□ 



dS 2 ) \0.2S 2 ) ^ \0AS 2 °-3St) \dW? 

where Wf,W 2 are independent Brownian motions. Here each asset is driven by two sources 
of noise. It may be convenient to rewrite the dynamics in a decoupled fashion: 

dS] = 0.3SI dt + axS} dW} 
dS 2 t = 0.2S t 2 dt + a 2 S 2 t dWf 

where Wf,W 2 are correlated Brownian motions. This may be simpler, because each asset is 
now driven by only one source of noise. 

The two things that we need to know are: 

(i) What are the volatilities 0"i,0"2? 

(ii) What is the correlation p between and W 2 ? 
Clearly, we must have 

<7 1 dW} = 0.1 dW] +0.4 dWl 
a 2 dW? = 0.4 dW} + 0.3 dWl 
Looking at the covariance processes, we must have 

a\ dt = (0.1 2 + 0.4 2 ) dt 
a\ dt = (0.4 2 + 0.3 2 ) dt 
G X G 2 p dt = (0.1 x 0.4 + 0.4 x 0.3) dt 

which are three equations in 3 unknowns, easily solved for <5"i,<T2,p: 

<7i = ||(0.1,0.4)|| 
(T 2 = ||(0.4,0.3)|| 

(0.1,0.4) • (0.4,0.3) 
P ~ ||(0.1,0.4)||. ||(0.4, 0.3)|| 

Note that the vectors on the right can all be read off the volatility matrix. 
Thus 

(0.1,0.4) -(W t \W?) 
||(0.1,0.4)|| 
~ 2 _ (0.4,0.3) -(W t \W t 2 ) 
||(0.4,0.3)|| 

It is clear that W 1 ,]^ 2 are continuous martingales. Moreover 

[W 1 ^ = t = [W% 

so that W 1 ,]^ 2 are indeed Brownian motions (by Levy's Characterization). Furthermore, 

[W\W 2 ] t = pdt 

as expected. 
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□ 



The above example can be generalized: 
Proposition 9.3 Give a multidimensional SDE dX t = b(t, X t ) dt + a(t, X t ) dWt, i.e. 

'dW}\ 





fb\t,X t )\ 




fo-u(t,X t ) . 


■ o- lrn (t,X t )\ 






dt + 






\dXfj 


\b n (t,X)j 




\a nl (t,X t ) . 


■ o- nm (t,X t ) J 



where Wt = (W/, . . . , W m )t is a standard m -dimensional Brownian motion. Let ai be the i 
row of the matrix a . Define 



tli 



W; = — — — for i = l,..., n 



Then (by Levy's Characterization) the W\ are n correlated Brownian motions, with correlation 

0~i ■ <Jj 



ll°i|| " ll^jll 

and we have dynamics 

dXi = b i (t,X t ) dt+\\o-i(t,X t )\\ dW\ fori = l,...,n 
Here each X 1 is driven by only one source of noise. 

Thus the "volatility" of a one-dimensional process of the form 

dX t = t idt + a 1 dW 1 + ... + a n dW n 



□ 



is 



U= ||(<7l,...,0|| = yfi* + • • • + a* 



What happens to the Black-Scholes PDE when we have correlated Brownian motions? 
Recall that this is 

00 + y rS n ^- + I T(aa tr )S n S m 1 -f^— - rC = 
f),qn ^ l > dS n dS m 



dt 



Note that (aa ) nm is just a n ■ a m , the inner product of the n and m rows of a. We have 
seen that 

Pnm 



0~n ' 0~m 



0~n ' 0~m 



0~n 0~m 



OnOr, 



and thus we obtain 

dC 
~dt 



+ E + \ £ PnrnO-nVrnS-S™^^ ~ ^ = 



where a n is the volatility of S* 
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10 Exercises 

1. An endowment option X is a very long term European call option. Typically, 

• At issue, the initial strike K$ is set to approximately 50% of the current stock price. 

• The options are inflation and dividend protected: 

— The strike price increases at the short term riskless rate. 

— The strike price is decreased by the size of the dividend each time a dividend is 
paid. 

• The payoff at expiry T is X T = (S T - K T ) + . 

We will make the simplifying assumption that the stock pays no dividends. This can 
be accomplished by regarding the stock price as the theoretical price of a mutual fund 
which starts off at one share, and reinvests all dividends in that share. We have, in the 
risk-neutral world, 

dS t = r t S t dt + a t S t dW t dA t = r t A t dt 

where S is the share and A the money market account (with Aq = 1). Clearly K t = K^Af. 
By changing the numeraire to A t , show that, when the volatility at is deterministic, 

X = S N(d + ) - K N(d-) 

where 

In ^ ± -a 2 T i f T 

d ± = — — TfjT— and a lv = ^ / of dt 

VavVT J- Jo 

2. Use the change-of-numeraire technique to show how to calculate the value of an option 
which pays the minimum of two assets S 1 , S 2 . Assume that the "real world" dynamics of 
the assets are Ito diffusions of the form 

dS\ = S l t [fii dt + di dwi] 

where m,(Ti are constants, and that the correlation of returns is a constant p. Further 
assume that SI has a continuously paid dividend with constant dividend yield q l . 

3. Consider a European call C on share S traded on FTSE. St and C are priced in pounds, 
but the strike of the call is in dollars. Initially, the option is at-the-money. The dollar 
strike does not change, but because exchange rates are not fixed, the pound strike does. 
Let X t be the ^fa-rate, Y t the 1 j^-rate. Assume dynamics 

dS t = a s S t dt + 5 s S t dWf 
dX t = a x X t dt + 5 x X t dW t x 
dY t = a Y Y t dt + 5 Y Y t dWj 

where W S ,W X ,W Y are correlated Brownian motions. 

3.1 Apply Ito's formula to show 

dY t = a Y Y t dt + 5 Y Y t {-dW t x ) a Y = -a x + 6%, 5 Y = 5 X 
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3.2 Let p be the correlation between W s and W x . Let W t = (W^W?) be a two- 
dimensional standard Brownian motion, and rewrite the above dynamics 

dS t = a s S t dt + Stas dW t 
dX t = OL X X t dt + X t ox dW t 
dY t = a Y Y t dt + Y t o Y dW t 

Show that we must have 

II— I|2 jr2 i i ||2 e2 1 1 ||2 c2 
ll^-Xll = °X ll°V|| = °Y \\ a S\\ = S 

&x ■ os = pdx^s &y ■ &s = -pSx^s 

3.3 The initial pound strike is Kq = So, and the initial dollar strike is K d = SqXq (at-the- 
money), which remains fixed. At maturity, the pound strike is Kt = K d Yr- Define 
Sf = S t X t to be the dollar price of S at time t. Show that 

dSf = S?[a s + a x + a s - ax] dt + Sf(a s + <r x ) dW t 

3.4 Now convert this to a system with a one-dimensional Brownian motion Vt- 

dSf = Sf[a s + a x +o-s- <t x } dt + S?5 sd dV t 

where 

S 2 sd = \\<rx + vs\\ 2 = (4 + $s + 2pSx$s) 

3.5 Now we have a plain vanilla call on an asset S d with (fixed) strike K d . Find the dollar 
price Cf of this option: 

C d = S d N{d + ) - e- rd{T - t] K d N{d_) 
where where is the riskless dollar rate, and 

ln§ + (r d ±±5 2 sd )(T-t) 



5 S d\/T — t 

3.6 Conclude that the pound price of the option is 

C t = S t N{d + )-e- r ^ T -^^N(d-) d± = SoXo , 2 X S 

\5\ + 5l + 2p5 x 5s){T-t) 



3.7 If we had tried to price the option directly in pounds, we would have had (explain 
this) 

C T = {S T - S (Y T /Y )) + 

Very naturally, we would have considered the numeraire Y t . This would have been a 
mistake, for although Y t is a traded asset (namely the pound price of a dollar note), 
this is not a non-dividend paying asset: Yt has a continuous dividend yield equal to the 
riskless dollar rate r^. Thus discounted Y t is not a Q-martingale. Instead, therefore, 
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consider the process Y t = Y t e Tdt . (i.e. all dividends = interest reinvested in the dollar 
money market account). Show that 

C t = Y t ^[{S t -K') + \T t ] 

where K 1 = e~ TdT Sq/Yq and Q is the equivalent martingale measure associated with 
Y. 

3.8 Find the Q-dynamics of St (with a two-dimensional standard Brownian motion). 

3.9 Convert this to S^-dynamics with a one-dimensional Brownian motion. 
3.10 Hence show that 

C t Y t [S t N(d+)-K'N(d-)] 

where 



d± 



ln§f 7 ±f6 2 (T-t) 



and 5 = 1 1 as — <Jy 1 1 = \J&$ + + ^pSs&y- 

3.11 Finally show that this coincides with the formula obtained earlier. 

In this case, you see that it is slightly easier to value the option in dollars than it is in 
pounds. 

4. Suppose the bond price dynamics are given by 

dp(t, T) = p(t, T)M(t, T) dt + p(t, T)v(t, T) dW t 
Show that in that case the forward rate dynamics are given by 

df(t, T) = a(t, T) dt + a(t, T) dW t 

where 

a(t, T) = v T (t, T)v(t, T) - m T (t, T) a(t, T) = -v T (t, T) 

[Hint: Apply Ito's formula to lnp(t,T), write this in integrated form, and differentiate 
with respect to T] 

5. Let {y(0,T) : T > 0} denote the zero-coupon yield curve at t = 0. Assume that, apart 
from the zero coupon bonds, we also have exactly one fixed coupon bond for every maturity 
T. enote the yield-to-maturity of the fixed coupon bond by ?/m(0,T). We now have 3 
curves to consider, the forward rate curve /(0,T), the zero yield curve y(0,T) and the 
coupon yield curve ?/m(0,T). 

5.1 Show that /(0, T) = y(0, T) + r ^ T) 

5.2 Assume that the zero yield curve is an increasing function of T. Show that in that 
case 

!/Af(o,r)<y(o,r)</(o,r) 

for all T. Show that the inequalities are reversed if the zero yield curve is decreasing. 
Explain this phenomenon in terms of simple economics. 
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5.3 Yield curves can be both upward and downward sloping. Can this be true for bond 
price curves p(0, T)? 

6. In the Cox-Ingersoll-Ross model, the risk-neutral short rate dynamics assumed are 

drt = (b — art) dt + a^/n dWt, a, b, a, ro > 



6.1 Explain (heuristically) why this process is mean-reverting and non-negative. 

6.2 This is an affine short rate model. By plugging p(t, T) = e j4 (*> T )-- B (*>' r ) r t ' m i Q the term 
structure PDE, show that we obtain two coupled ODE's 

B t = aB+^a 2 B 2 -1 B(T,T) = 
A t = bB A(T, T) = 

6.3 To solve the Riccati equation for B, try a solution of the form 

X(t) 



B(t,T) 



cX(t) + d 



Choose c to ensure that a + \o 2 — c 2 = 0. Show that we then obtain a order linear 
differential equation 

X t + kX = —d where k = — a + 2c = \f a 2 + 2a 2 

6.4 Solve the ODE to obtain 

X(t) = -[e^ T -^ -1] 

K 

6.5 Hence show that 

2 ( e K(T-t) _ d 

B ^ T ) = o rA v nvr-t) u where K = v 7 " 2 + 2^ 2 

2k + (a + Kj(e K( --' ^ — 1) 



6.6 Verify by differentiation that 



2b 

A(t, T) = — 2 In 



2 Ke |(«+«)(T-t) 



2 K + (a + K )(e K (' r -')-l) 



7. 7.1 Show that the Hull- White model dr = (6(t) — ar) dt + a dWt is obtained if one starts 
with a HJM model given by 

df(t, T) = a(t, T) dt + ae" a(T -* ) dW t 

Hence compute the function 9(t) which will make the short rate model fit the initial 
term structure: 

9(t) = m0,t) + af(0,t) + ^-[l-e- 2at } 

where |/*(0,T) : T > 0} is the observed term structure of forward rates. It follows 
that the Hull- White model can also be fitted to any initial term structure. What is 
the distribution of the forward rate f(t,T)? 
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7.2 Show that bond prices in the Hull- White model, fitted to the initial term structure, 
are given by 



P(t, T) = exp (/(0, t)B(t, T) - -B 2 (t, T)(l - e~ M ) - B(t, T)r t 

where B(t,T) = £[1 - e~< T -% 

[Hint: The Hull- White model is an affine term structure model, i.e. p(t, T) = e A ^ t ' T ^~ B( - t ' T ^ Tt . 
B(t, T) is readily calculated. We can now find 

A(t,T)= [ -6(u)B(u,T) du+\o 2 B 2 (<d,t) 
Jt 2 

where 9(t) is as in (a), i.e. b(u) = e~ au ^x(u)e au , where x(t) = e~ at ro + f* e ~ a (*-") du. 
Integrating by parts leads to 

A(t,T) = f(0,t)B(t,T)+ln P{1) i] 



p(0,t) 



+ 



a 2 



J b 2 (u, T) - B 2 {0, u) du + B 2 {0, t)B(t, T) 



which simplifies to give the required result.] 

7.3 Show that the Hull- White price of a call C with strike K and maturity T on a bond 
p(0,S) (where S > T) is given by 

p(0,S)N(d + )-Kp(0,T)N(d-) 

where 

, 111 Kp(0,T) ^ 2 a av 1 

d± = 



and where 

fT J2 



a 2 av T = f ^(e- aS - e' aT ) 2 e 2at dt = ^(1 - e -< s - T )) 2 {l - e~ 2aT ) 
Jo o 2a 

8. Consider the domestic and the foreign bond market, with bond prices denoted by p d (t,T) 
and pf(t,T) respectively. Take as given a standard HJM model for the domestic forward 
rates f d (t,T) 

df d (t, T) = a d (t, T) dt + a d (t, T) dW t 

where W% is a multidimensional Brownian motion under the domestic martingale measure 
Q. The foreign forward rates are denoted by ff(t,T), and their Q-dynamics are given by 

df f (t,T) = a f (t,T)dt + o- f (t,T) dW t 

Note that the same Brownian motion drives both bond markets. The exchange rate X (in 
units of domestic currency per unit of foreign currency) has Q-dynamics 

dX t = n(t)X(t) dt + a x (t)X(t) dW t 
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Show that under the domestic martingale measure the foreign forward rates satisfy the 
modified HJM drift condition 



a f (t,T) = a f (t,T) 



af(t,s) 



ds - G x (t) 



9. A common implementation of the HJM framework uses the following forward rate dynam- 
ics: 

df{t,T)=a(t,T) dt + (a 1 ,a 2 e-^ T - t '>) • (dW 1 (t),dW 2 {t)) 

where oi, a 2 , A are non-negative constants, W\, W 2 are independent Q-Brownian motions, 
and Q is the equivalent risk-neutral measure. 

This is a two-factor model. The first factor W\(t) can be interpreted as a source of noise 
that lasts a long time, affecting all maturities equally. The second factor W 2 {t) affects 
short maturity forward rates more than the long term rates (why?), and thus adds some 
extra volatility to the short term rates. 

9.1 Show that the HJM drift conditions imply that 

a(t,T) = af(T - t) - ^e"f (T-t) {e -^T-t) _ 1} 

A 

9.2 Hence show that 

f(t, T) =/(0, T) + o\t{T - t/2) - 2(a 2 /X) 2 [e~ XT (e xt - 1) - 2e -( A /2)T (e (A/2)t _ ^ 

+ a 1 W 1 (t)+a 2 [ e ~(V2)(T-«) dw ^ 
Jo 

9.3 Show that the spot rate follows the process 

rit) = f(0, t) + \a\e - 2(a 2 /A) 2 [l - e~ <*/*>']=> 

+ axW^t) + a 2 e' {x ^ 1 f e^ u dW 2 {u) 

Jo 

9.4 Is the short rate a Markov process, a Gaussian process, a stationary process? Explain. 

9.5 Calculate the price C{t) of a call option on the zero coupon bond p(t, T). Assume that 
the option has strike K and expiry r, where t < r < T. 

[Hint: Let p(t,r) be the numeraire. You know the HJM dynamics of zero coupon 
bonds under Q, so the dynamics of p{t,T)/p{t,r) under the EMM for p(t,r) should 
be easy to find. Of course, something is going to be lognormal. Now use the general 
option pricing formula.] 

9.6 As a check, assume that o\ = 0.2, a 2 = 0.3 and A = 2. Calculate the value of a call 
option on a two-year zero coupon bond with strike 0.9 and expiry 1 year. Today's 
prices are P(0, 1) = 0.9,P(0,2) = 0.81. I get 0.076 (but I could be wrong, of course). 

10. Consider a convertible bond X which, at To, allows the owner to convert the bond to c 
shares S of common stock. The bond is a zero coupon bond with face value 1.00 and 
maturity T\ > Tq. The aim of this problem is to find the price the convertible bond at 
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some future date t < Tq. We will model the short rate using Ho-Lee dynamics. Initially, 
the (instantaneous) forward rate curve is flat with /(0, T) = tq for all maturities T. 

We work under a risk-neutral measure Q where the share has dynamics 

dS t = r{t)S t dt + a s S t dW t 
and the short rate has dynamics 

dr(t) = d(t) dt + o r dW t 
Here Wt is a two-dimensional Q-Brownian motion, and crs,cr r are constant vectors. 

10.1 Let p(t, T) be a non-convertible zero coupon bond with face value 1.00 and maturity 
T years. Calculate the observed term structure of bond prices {p*(0,T) : T > 0}. 

10.2 Let Q be the forward risk-neutral measure for maturity T\ years (i.e. the EMM for 
numeraire p(t,Ti)). By decomposing the convertible bond into its option and bond 
parts, show that 

X = cp(0,Ti)E [(S To - + P (0,T 1 ) 

where & = flffa- 

10.3 The Ho-Lee model is an affine term structure model, i.e. bond prices are of the form 

p(t,T) = e Mt,T)-B(t,T)r(t) 

By substituting this expression into the term structure PDE, show that 

rT 



B(t,T) = T-t A(t,T)= [ 6{u){u-T) du+\a 2 r .{T -tf 

Jt 6 



10.4 Show that 

d ^~ = r(t) dt-a r (T-t) dW t 

10.5 Fit the Ho-Lee model to the initial term structure of forward rates: Show that 

0(t) = CT 2 r t 

10.6 Hence show that 

r(t) =r + -a 2 t 2 + a r W t 

10.7 Hence, using the known initial bond prices and short rate dynamics, show that future 
bond prices are given by 

p(t,T) = e -k a rt(T-t) 2 -r(t)(T-t) 

What is the distribution of p(t, T) under Q? 

10.8 Show that St has dynamics 

^ = (a s + a r (T 1 -t))dW t 

under the measure Q, where Wt is a Q-Brownian motion. 
What is the distribution of St under Q? 



Interest Rate Modelling 



89 



10.9 Deduce that 
where 

d± = 



X = cS N(d+) -p(0,Ti)JV(d_)+p(0,Ti) 



1 r Tl 



\a s + <Tr(7i - s)|| 2 ds 



11. The aim of this problem is to calculate the price of an in-arrears caplet in the Ho-Lee 
model, where the short rate has riskneutral dynamics 

dr t = 6{t)dt + o dW t 

Here, a is a constant, and Wt is a 1-dimensional Brownian motion under the risk neutral 
measure Q. The caplet has payoff 

0.5max{L- i? c , 0} 

at expiry = 1 year, where L is the 6-month spot LIBOR rate in 6 months' time, and R c 
is the cap rate. Use the following data: 



P(0,T) 


e -r T 


^0 


10% 


R c 


12% 


a 


10% 



Here P(0, T) is the default-free zero coupon bond with face value 1 and maturity T. 



We proceed as follows: We first show that the caplet is equivalent to a portfolio of put 
options on zero coupon bonds. Then we recast the Ho-Lee model within the HJM frame- 
work in order to fit it to the observed (flat) term structure, and calculate the prices of zero 
coupon bonds. Finally, we calculate the prices of vanilla options on zero coupon bonds. 

11.1 First show that a caplet can be regarded as a portfolio of 6-month put options on 
the 1-year zero: 

Caplet =(1 + R C AT) put options on P{t,T-z) with strike 



1 + R C AT 
and maturity T\ 

where T x = 0.5, AT = 0.5, and T 2 = T x + AT = 1. 

11.2 The Ho-Lee model is an affine short rate model, with bond prices of the form 
P(t, T) = e Mt^)-B(t,T)r t _ g y substituting this form of P(t, T) into the term structure 
PDE, show that 



B(t,T) = T-t 

A(t, T) = - f 9(u)(T -u)du + \a 2 {T - tf 
Jt 6 
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11.3 In order to fit the short rate model to the observed term structure, we recast it in 
the HJM framework. You may use the facts about the HJM model which are stated 
on the formula sheet. 

Using the relation between forward rates and zero coupon bond prices and the value 
of B(t, T), show that the instantaneous forward rate f(t, T) has a constant "volatility" 
a, i.e. that the forward rate dynamics are 

df(t, T) = a{t, T)dt + a dW t 

for some function a(t,T). 

11.4 Use the HJM drift conditions to show that a(t,T) = a 2 {T - t) 

11.5 Hence show that 

r t = r + \a 2 t 2 + aW t 

and conclude that 

dr t = o 2 t dt + a dW t 

11.6 Next, show that 

A{t,T) = - l -a 2 t{T-t) 2 
and thus that zero coupon bond prices are given by 

11.7 Now that bond prices have been found, we will price bond options. Recall the general 
option formula stated in the Formula Sheet. Change the numeraire to the 6-month 
zero coupon bond P(0, Ti). Let P t = p^'^j • Write down the dynamics of Pt under 
T\ -forward neutral measure Qi. 

11.8 Show that Pr 1 is lognormally distributed under Qi, and find its distribution. 

11.9 Similarly, find the dynamics of Pt = i under the T2-forward measure Q2. Show that 
Pt 2 is lognormally distributed under Q2. 

11.10 Deduce the following formula for a call option on P(t, T2) with strike K and maturity 
Ti: 

C = P(0,T 2 )N(d + ) - KP(0,T 1 )N(d^) 

Write down expressions for d±. 

11.11 Use put-call parity and the table of the normal distribution to find the price of the 
original caplet. 



